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FOREWORD 


This report vas prepared hy Lockheed Missiles and Space Conipany under NASA 
Contract NASW-2553* The period of performance of the research reported was 
September 1973 throu^ April 197^- Technical direction for this contract 
originated from The Director, Physics and Astronomy (Code SO), WASA Head- 
quajTters^ H. D. Calahan, Contract Monitor. 

•nie contract was managed by Robert Weiss, Guidance and Control Systems 
organization of the I/tSC Space Systems Division. Chief investigator and 
principal contributor to the report was A. D. Anderson, Atmospheric Physics 
Branch, Radiation Physics Laboratory, of the IMSC Research/Development 
Division. 

This Final Report, LMSC Report d 4-2093^ (Guidance and Control Systems report 
GCS/ 39 ^/^ 11 ) docximents and summarizes the results of the entire contract 
work, including concliisions and recommendations based upon the results 
obtained. It incorporates the material presented in the two Quarterly Progress 
Reports previously submitted: 

QPR No, 1 IMSC/D38lt420, GCS/3475/6211, 73 Dec 12 
OPR Ho. 2 LMSC/D 387326 , GCS/ 3529 / 6211 , 7*^ Mar 12 
This final report is divided into three parts. Part I, Definition of the Problem, 
cont iins a detailed description of the problem and of the Solar and Geophysical 
Indices used in the investigation. Part XI, Analysis, presents an analysis of 
five stratospheric warming cases. Fart III presents an hypothesis concerning the 
relatlonsld.p between auroral electrojets and stratospheric warmings. 

The authors wish to thank Drs. John Evans and Robert Vemey of the 124SC Atmospheric 
Physics Branch for their helpful suggestions and guidance during the study. 



StMHART 

The primary objectives of the study are to establish the details of 
the stratospheric warming processes as to time, area, and intensity and 
to. try to correlate the warmings with other terrestrial and solar phenomena 
occurring at satellite platform altitudes or observable from satellite plat- 
forms. The main purpose of the study is to find if there are links between 
the perturbed upper atmosphere (mesosphere and thermosphere) and the strato- 
sphere that could explain stratospheric warmings. 

Stratospheric warmings are large-area phenomena containing great amounts 
of energy. They probably have major influence on the weather over large 
areas, on observations from satellites and high -altitude platforms, on the 
dynamics of the upper and lower atmospheres, on atmospheric chemistry, on air 
pollution observations and on backgrounds for surveillance and earth resources 
observations. Little correlative data has been brought to bear on defining 
the relationship between energetic events above 50 km and the warmings at 
about 40 km. Current rocket and balloon-borne IR sounding instruments do 
not provide good coverage on a global basis. It is difficult to provide 
dynamic readings throughout a warming by means of ground based instrumenta- 
tion. Satellite and other high altitude data, therefore, need to be applied 
to the problem. 

Although meteorological analyses indicate that the warmings result from 
the transfer of energy from the troposphere to the stratosphere, the 
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initiating or trigger mechanism is unknown. In other words, no one has 
attempted to predict stratospheric warmings. The emphasis of this effort 
is to attempt to find clues to possible trigger mechanisms arising in the 
upper atmosphere or solar activity. Several steps are Involved: (a) Deter- 

mine from the literature possible correlations of warmings with geophysical 
events; (b) Find solar and geophysical indices that represent the pertinent 
events; (c) Define the stratospheric warming events (start, raEiximum, end) 
so that these Indices may be applied to uncover the trigger mechanism(s). 
Existing data and literature have been searched to find the identified cases 
of stratospheric warmings. Case histories of these events have been chosen 
by assembling correlative high altitude satellite, rocket and ground based 
measurements, geophysical characteristics and solar activity data. Existing 
data from NIMBUS spacecraft sensors such as the MRIR, SIRS, IRIS, and SCR ' 
were used to define the existence and characteristics of the warming in the 
ao-50 tan-altltude regime. Data* from various particle and radiation experi- 
ments on satellites were used to determine some of the energetic conditions 
existing in the upper atmosphere prior to and during the warming events. 

The study of six different stratospheric warmings has Indicated that the 
magnetic indices A^, AE and erfiibit a marked increase 1 to 3 days 
preceding the onset of- a warming at 10^^ (30 km), and a marked decrease on 
the day of onset and k days preceding. It Is proposed that eleetrojet 
heating forms a "cidtical level" In the upper mesosphere, and that planetary 
waves propagating from the troppsphere Interact with the critical level, 
releasing heat at that level. Waves reflected by the critical level release 
heat^ in the stratosphere by interference with upward-propagating waves seen 
as a warming. 
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The stiidy Is divided into three parts* In Part I, stratospheric warm- 
ings are described, the correlation problem is . defied, and the various 
geophysical indices selected for the analysis are discussed. In addition, a 
bibliography on stratospheric warmings in included. In Part II, case his- 
tories for recent stratospheric warmings are given. The day-to-day detailed 
analyses of the warmings and geophysical Indices are described. Conclusions 
are presented concerning the correlation between the geophysical indices and 
the onset of stratospheric warmings. Part III presents a brief hypothesis 
concerning the relationship between auroral electro Jets and stratospheric 
t^arml2}gs . 
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1. STRATOaPHERIC WARMINGS 

i ■ 

k c^namic event in the stratospheric circulation which is princi- 
pally characterized hy a temperature increase in higher latitudes immediately 
above the stratonull level greater than over a period of ten days 

or less, accompanied by a disruption of the usual westerly zonal circum- 
polar flow of the stratospheric winter circulation. Is called a Strato- 
spheric Warming* 

Stations centered under the warming region as it sweeps poleward 
provide a very dramatic record of the magnitude of events idiich must be 
in progress in the upper atmosphere. Temperature increases of from 30C* 
to 90c “ over a period of a week are the rule, and stratospheric wind 
changes of the order of 100 meters/sec out of the west to 50 meters/sec 
from the east are comaon (Webb, 1966). 

Figure 1 contains a sequence of rocket soundings from WestGelri22ish,Gernany, 
taken during a major warming In the winter of 1967-68 (Labltzke,l97l) . Profile 1 
indicates the beginning phase of the warming with the maximum waxming at 
60 km (to +11*^3) In Profile 2, the maximum wanning has descended to 46 km and 
Increased to tl8^* Profile 3 (21 Dec. I967) the maximum warming is at 
43 km and the temperature is +40*tJ, Profile 4 represents a decreasing 
phase; the maximum wanning at 44 km has decreased to -4*’C. Profile 5 is an 
avereige atmospheric temperature profile from Groves (I97O). Comparison 
of Profile 5 with the other profiles indicates that the warming affected 
a depth of atmosphere of approximately 34 km. 
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Stratospheric wannlngs are large-area phenomena containing large 

amounts of energy. They probably have major Influence on the weather 

over large areas, on observations from satellites and high-altitude plat- 

* 

forms, on the dyna m ics of the upper and lower atmospheres, on atmospheric 
chemistry, on air pollution observations and on backgrounds for surveillance 
and earth resources observations. 

Little correlative data has been brou^t to bear on defining the 
relationship between energetic events above 50 Ion and the warmings at 
about 40 km. Current rocket and balloon-borne IR sounding instruments 
do not provide good coverage on a global basis. It is difficult to pror 
vide dynamic readings throughout a warming by means, of ground based 
instrumentation. Satellite and other high attitude data, therefore, need 
to be applied to the problem. 

Although stratospheric warmings have been studied extensively since 
their recognition in 1952, their characteristics are not well known or 
established. They occur mostly in winter and early spring seasons and 
show largest temperature changes in the northern and southern high-latitude 
regions. They start at 50 km or above and later sometimes appear at lower 
levels, penetrating to 25-35 km. The energy density of direct solar 
radiation on an element of atmosphere is smaller by many orders of magni- 
tude than that req^ulred to explain the atmospheric temperature rise during 
a warming. Integration of this energy over time and area would be required 
to explain a warming in terms of direct heating. The heat input into a 
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square cm column of the stratosphere "between 30 and 50 km altitude needed 

o' 9 ■ -2 

to raise the temperature 50 K amounts to some 5 x 10 ergs cm . Reported 

15 2 

warmings may cover 10 cm so the total energy involved can be in excess 

24 23 

of 10 ergs. This energy might be compared with 5 x 10 ergs reported 

in the solar storm of 11 February 1958 in the red arc aurora, and a cora- 
puted total solar particle energy of 10 ergs in such a storm. Muench (1965) 
considered the energy transformations during the stratospheric warming of 
22-24. January 1958» His analysis indicated that the energy increase was 
caused almost entirely by a transfer of energy from the troposphere to 
the stratosphere. A report by Hirota (196 t) on an observed upward propa- 
gation of stratospheric warming events supports this conclusion. Matsuno 

(1971) treated the interaction of vertically propagating planetaiy waves 
from the troposphere with the zonal flow in the stratosphere. His model 
produces many observed features present In stratospheric warmings. 

Trenberth (1972) used a model of this type to study stratospheric warm- 
ings. His model reproduced realistic stratospheric warmings in late 
winter. Although Trenberth neglected the chemistry, Clark (1970) had 
previousiy used a similar model (with somewhat fewer spherical harmonics 
and coarser vertical resolution) which did Include the simple Chapman 
scheme of oxygen photochemistry. Clark was able to reproduce numerically 
the observed tendency for the large-scale motion field to shift, the maxi- 
mum of vertically integrated ozone poleward from the low-altitude source 
region. In such a model, the chemistry interacts with the dynamics via 
the heating rate Q, which will depend not only on infrared radiation and 
a parameterization of small-scale convection in the troposphere, but also 
the absorption of radiation by the variable amount of ozone. 
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A survey of recorded ina^or and minor stratospheric -warmings shows 
that in January of every year between I966 and I971 at least one event 
occurred In the Northern Hemisphere between 45 ^ and 80 ° latitude. These 
same years produced events in the Southern Hemisphere in July and August. 



2. THE PROBUM OP DEPINIHG A GIVER STRATOSPHERIC WARMING EVENT 


A given stratospheric warMng (SW) can he defined provided that the 
date and time of its start, maximum and end can he determined from the 
temperature at all aJ-titudes that it affects. One major prohleni with char- 
acterization of Slfe occurring in the past is that they occur at altitudes 
between about 20-65 kni (Fig. l). Until recently, the higher altitudes 
were only accessible to measurement with meteorological rockets. There are 
several disadvantages for defining warmings with rocket data. Although 
there are currently more than 20 meteorological rocket launch sites, these 
are. mainly land-based in the Northern Hemisphere. Hence, the coverage is 
uneven. Also, the launches are made closer to a weekly than daily basis. 
However, even if daily launchings were made, this would be insufficient to 
characterize warmings, Inasmuch as the warming centers probably would not 
lie over the rocket launch sites. In addition, the warming centers move . 

The coverage and launch frequency of radiosondes carried by balloons 
is much better. There are at least 20 times as many radiosonde launch 
sites than rocket sites. The coverage is best In the Northern Hemisphere 
over land; the coverage is much poorer in the Southern Hemisphere because 
of the greater amount of water area. Also, the launch frequency (one to 
four times daily) is much greater than rockets. The main disadvantage is 
that only the lower altitudes at which SWs occur are accessible to the 
balloons (ceiling altitude about 30 km). The 10-mb surface (3l km) is the 
hipest pressure altitude for which hemispheric charts, derived mainly from 
radiosonde data, are available from the National Weather Service on a dally 
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basis for the Northern Hemisphere warmings occurring during December I969 - 
January I97O, December 197^ ” January 197Jj and January 1973 " February 

1973. 

Ihe real hope for the future lies in satellite remote sensing tech- 
niques that scan worldwide twice daily. One example is the NIMBUS SIRS 
instrument data from which temperature profiles can be derived for the 
altitudes of interest on a worldwide hasis. At present, there are some 
problems in obtaining data at the hi^er altitudes above km. Therefore, 
at this time,, the analysis of a given SN must be based on radiosonde, rocket- 
sonde, and satellite (radiance) data* The smallest time interval that can 
be used in characterization is about one day. It is not possible to specify 
warming start, maximum, pr end at a given altitude in a time interval less 
than about one day. The, definition must be Incomplete for the warmings 
considered in this report because little data, is available at the higher 
altitudes above 40-50 km. The radiance measurements made by the SCR on 
Nimbus 4 provide radiation information from the highest levels in the 
atmosphere so far obtained. Descriptions of the instrument can be found 
in the literature (Houghton and Smith, 1970) and the method used by the 
Oxford University Group to derive temperature profiles from the six radiances 
using a maximum probability technique has been described by Rodgers (1970)- 
A aeries of daily charts for levels up to 1 mb (48 km) were constructed 
to cover the major northern hemisphere warming period in the middle of 
February 1972. 
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3* POSSIBLE CORRECiAaPlON OP WAHMISGS WITH GEOPHYSICAL EVENTS 

S'bra'tospheric winnings favor the winter awd early spring seasons; 

they exhibit the greatest change of temperature In the polar regions. 

They start at 50 km or above and later sometimes appear at lower levels. 

In some cases, the warmings do not penetrate downward as far as 30 km. 

* 

These very high level relatively limited stratospheric warruirigs may roipre- 
sent the effects of comparatively limited penetration of the atmosphere 
by solar corpuscular emanations of only light to moderate intensity, 
which should be expected to occur not infrequently. 

The major warmings of February ' 1952 and January 1963 , which came 
well after the years of hipest sunspot activity, both followed sudden 
bursts of strong geomagnetic disturbance with a moderate increase' of sun- 
spots, and In 1952 ,wlth a strong auroral display. The geoinagnetically 
quiet years of 1953-1955 and I964-I966 were apparently strikingly' lacking 
in stratospheric warmings. The fact that of the sixteen years, 1951 to 1966 
inclusive, the only six years lacking stratospheric warmings were also 
the six geomagnetically most quiet years would seem significant. With the 
exception of one weak warming at the beginning of March I956, all 14 
warmings, plus the two occurring in February I952, originated during the 
months of Januaiy and February, when Northern Hemispheric auroral activity 
is very high, rising to a presumed early March peak frequency which may 
be even earlier in arctic latitudes. It is to be noted that none of them 
occurred during December, the month of the presumed winter solstice minimum 
of auroral frequency, and definitely the ccOdefft month in the arctic strato- 
sphere. During the period I95I-I966, all of the observed sudden 



stratospheric vannlngs developed on the northern geomagnetic pole side of 
the geographic pole, i,e., essentially in the equatorward dip or trough 
of the zone of maxiraum auroriO. frequency. 

The correlation of the time and space distribution of auroral activity 
vith sudden changes of stratospheric tenq)erature indicate a possible 
extraterrestrial control, perhaps solar corpuscular invasions of the upper 
atmosphere i^ich produce auroras and strong geomagnetic disturbances. If 
solar corpuscular invasions are to cause and control the entire syndrome 
of associated auroral, and thermal phenomena, then they should possess the 
follovlng properties and characteristics: 

1) The energy contained in one Of the strongest corpuscular Invasions 
must be adequate, at whatever level of efficiency its conversion Into ther- 
mal energy takes place, effectively to trigger the change of the temperature of the 
entire polar stratosphere in a relatively short time and in the virtual absence 

of sunlight from normal midwinter to normal midsummer conditions, such as 
occurred in the massive stratospheric warmings of January 19^3. Less intense 
solar events (corpuscular Invasions) will be restricted in their significant 
thermal effects to above the 50 km level. 

2) There must exist sone fairly effective means of downward trans- 
mi^islott in the upper atmosphere of the direct thermal and magnetostrophic 
wind effects of strong solar corpuscular penetration of the higher atmo- 

f ' 

S ' 

sphere to at least the 50 -km level. 
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3) Strong solar corpuscular penetration of the upper atmosphere 
should be favored^ either as to frequency of occurrence or level of lovest 
penetration in the atmosphere, or both, on the geomsignetic pole side of 
the hemisphere In the general region where the zone of maximum auroral 
frequency extends to its lowest geographical latitude, 

4) The initial burst phase of a massive warming of the middle or 
upper stratosphere should be roughly comparable in areal extent to that 
of a major auroral display (strong solar corpuscular penetration of the 
upper atmosphere). 

Basically, it Is the first two of the required characteristics of 
solar corpuscular invasions as noted above that are essential to their ■ 
operation as a trigger . of sudden massive stratospheric warmings. 
They must bring a large part of the necessary energy with them, and this 
energy must be transmitted effectively downward from the level of Its direct 
absorption in the higher atmosphere 7 O-IOO km) to the middle or lower 
stratosphere. 

The forms of corpuscular radiation which merit attention for present 
purposes are: 

1) Electrons, in the energy range 1-100 keV, (and also in the range np 
to 1000 keV) which precipitate over the auroral zone. 

2 ) protons, in the energy range 1-1000 MeV, originating in solar 
flares, and precipitating over the geomagnetic polar cap. Included in this 
latter group may be a proportion of a-particles, and possibly heavier 
elements, 
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Protons in the energy range 1-100 keV normally do not penetrate below 
100 kmj and hence are not considered as main phenomenon. All these phenomena 
tend to occur within the higher geomagnetic latitudes hounded hy about 
magnetic latitude 45 ^ and mainly around 65 - 70 *^# but during major disturbances 
of the earth' s, magnetic fields extending a few degrees nearer the equator. 

The depth of penetration of electrons into the atmosphere can be 
evaluated from known characteristics of energy loss, and by adoption of 
a model atmosphere. For example, a 10 keV electron may penetrate to a 
minimum altitude of 98 km, while a 400 keV electron will reach 64 km. 

There remains the question as to the fraction of total energy which is 
available for heating the atmosphere, Chamberlain has estimated a fraction 
of about as appearing as ■ heat. Dalgarno estimates the fraction as about 
2C^, but does not exclude a value of 5 OS 6 . The heat input from electrons -In 

auroral latitudes and at altitudes of 90-120 km is thus of the order of 

-2 -1 

1 erg cm sec for solar maximum and for zonal and long terra average flux. 
After major solar flares, the corresponding zonal average is probably at 
least one order of magnitude greater for periods of a few days. . 

In contrast to the zonal character of auroral electron precipitation, 
the incidence of soleir protons is distributed over the polar caps. Incom- 
ing protons are assumed to move in Stomer trajectories, and the magnetic 
rigidity necessary to permit entry at a given angle to the vertical and 
at a given latitude may be calculated provided the spatial character of 
the earth's field is laiown. As a rough guide, the precipitation of solar 


IJ.0 



protons inac/ be considered as bounded by the auroral zones, but in major 
events, marked effects occur as low as 50 ° geomagnetic latitude. 

The lowest observed ionization from solar protons has been found at 
42 km in the southern polar region during winter darkness. Loss of energy 
from an incident proton is continuous along the trajectory, but is essen- 
tially concentrated over the last few kilometers of path. Reid has calcu- 
lated that a 50 MeV proton will lose of its energy between altitudes 
of 50 and 42 km, between 44.5 ahd 42 km, arid XOjh between 42.1 and 42.0 km. 
Hence, in a particular solar proton event, the increment of energy released 
in a given altitude range depends upon the spectrum of the flux and the 
geomagnetic latitude. An upper boundary of about 100 km is typical for 
most events, with release of energy down to a probable lower limit of 30 km 
in the largest events. 

The main effects of particle precipitation are those which have been 
used extensively to determine the morphology of the precipitation zones, 
viz., optical excitation and electron- ion pair production. The effects of 
interest are the heating produced directly by the particle flux, and the 
possibilities of indirect effects. Among the latter may be included the 
modification of minor constitutuent concentrations, and the initiation 
of atmospheric wave motion. The conversion of the kinetic energy of inci- 
dent particles into heat has been discussed by Bates and Chamberlain and 
the values of conversion efficiency used are between O.I 5 and 0,5. 



The major temperature increases are naturally associated with the 
largest fluxes of electrons in aurorae. Electrons in the energy range of 
1-10 keV affect the altitude range I5O-IOO km, and here the increases in 
temperature during more Intense fluxes occur at rates up to lOOK min"^. 

In general, due to the many factors which enter Into a quantitative dis- 
cussion of temperatures of' the higher atmosphere, and also due to the 
difficulties of observation of this matter, the effects on the mesosphere 
of particle heating are not well known. 

- The possibility that particle Influx _ml^t produce changes in con- 
centrations of minor atmospheric constituents, and that the return to 
equilibrium of these concentrations might be on a longer time scale than 
the duration of Intense fluxes, has been considered by a few workers. The 
dissociation of molecular oxygen by auroral electron and proton impact has 
been investigated by Maeda in relation to the warming of the polar meso- 
sphere. He found electrons with a differential spectral distribution 

-4 -8 - 2-1 

of the form E , or a little over 10 cm sec at 80-85 km altitude. 

Protons, distributed' according to a differential spectrum of E” , gave 
similar results, though larger values than electrons below 80 km. This 
rate coefficient is nearly the same as that due to ultraviolet dissociation 
of Og In the Schumann-Runge continuum (I40O-I750A) , and the height of maxi- 
mum dissociation is nearly the same, viz., 90 km. Thus, the process of . 
dissociation by particles may be significant during large fluxes, especially' 

I 

in polar winter darkness. 
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The generation of traveling atmospheric vaves is one consequence of 
particle precipitation, presumably due to a combination of spatial and 
temporal changes, including pulsations, In the incident flux, with corres- 
ponding changes in local heating. King refers to the evidence for these 
waves, which in some modes of propagation proceed from the auroral zones 
towards mid-latitudes with velocities ranging between 400 and TOO m sec“^. 
Other workers have considered the heating aspect of these waves. Thus, 
these disturbances, though clearly associated with the auroral zone Initially, 
may be included in the general subject of atmospheric wave motion, and the 

heating induced thereby. Especially Inqwrtaat Is the Jo\ole heating associated 
with the auroral electro Jets discussed later. 



k, SOLAR ARP GEOPHYSICAL IRPZCE3 


There is a wealth of operational data on solar and geophysical indices.. 

The main problem is deteiminlng which data are pertinent for correlation 
analysis involving the upper atmosphere and the warmings. In view of the 
fact that no direct correlations have been found between stratwarmings 
and the upper atmosphere phenomena represented by the indices^ the selection 
of the Indices to be used must be made on the basis of a possible connection 
and not a probable one. In other words, we are dealing with a strictly 
exploratory situation inasmuch as no connection may exist between strat- 
wannings and higher atmosphere phenomena. Althou^ meteorological research 
to date indicates that stratospheric warmings are caused by a transfer of 
energy from the troposphere to the stratosphere (Matsuno, 1971) ^ neverthe- 
less the possibility has not been completely ruled out that such events are 
triggered or enhanced by certain conditions occurring in the upper atmo- 
sphere above 80 km. 

The solar .and geophysical indices used are listed in Table 1; in-house 
values are indicated by a X mark for the months containing warmings. These 
particular indices were selected on a preliminary basis because they all 
appear to be correlated with upper-atmospheric heating phenomena. All 
of the indices except the Hell flux and D^_j_ are available from the Environmental 
Data Service, NOAA, Rational Geophysical and Solar-Terrestrial Data Center, 
Boulder, Colo., in the monthly publication "Solar-Geophysical Data. " Each 
index in Table I will now be discussed in turn. 
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PRICES qSED TOR STRATWftBM ANAiySIS 


INDEX 

Apr 

May 

June 

July 

Aug 

Dec 

Jan 

Dec 

Jan 

Jan 

Feb 


*69 

'69 

*69 

'69 

^69 

*69 

*70 

*70 

*71 

*73 

'73 

Solar Flares 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

Dally Flare Index 






X 

X 

X 

X 

X 

X 

S.I.D. 

X 

X 

X 

X 

X 

X 

X 

X 
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Solar Flares 
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Thfi power emitted by the sun is -- 4 x 10 ergs sec , of which 
24 -1 

^ 5 X 10 ergs sec , or 5 million horsepower per sq. mile, irradiates the 

atmosphere of the earth. Ionizing x-rays and ultraviolet radiation, 

which constitute only ^ lo"^ percent of the solar flux 5 x 10^^ ergs 

sec are absorbed in the outer atmosphere to produce the ionosphere. A 

20 -1 

considerably greater amount of energy, ~ 10 ergs sec , arrives at the 
magnetosphere in the form of solar-wind plasma, even under solar-quiet con- 
ditions. When the sun is intensely active, the energy delivered in ionizing 
radiations, both electromagnetic and corpuscular, can increase by orders 
of magnitude. No comparable energy supplies are available elsewhere in the 
sun-earth‘ system. 

32 

A great flare on the sun can release 10 ergs in x-ray ultraviolet 
and energetic-particle radiation. Every link in the system from sun to 
earth reacts to such a catastrophic event. After interacting with the inter- 
planetary medium, the magnetosphere, and the upper atmosphere, the last 
residue of energy from the sun’s explosive gift to earth is a faint red 
airglow that suffuses all the sky from north pole to south pole. Although 
the glow is composed only of the single, characteristic emission line of 
oxygen at 6300A, following a great flare the full sky may radiate 10 ergs 
of this monochromatic energy. 

A flare is a short-lived sudden increase in the intensity of radia- 
tion emitted in the neighborhood of sunspots. It is best seen in Hot, 
and usually occurs in the chromosphere. Flares are characterized by a 
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rise time of the order of minutes and a decay time of the order of tens of 
minutes. The total energy expended In a typical flare is about 10^^- ergs; 
the magnetic field Is extraordinarily high, reaching values of 10^ to 10^ 
gauss. Optical flares are usuaJ-ly acccajapanied by radio and x-rsy bursts, or 
flares, and occasionally by high-energy particle emissions. The optical 
brightness and size of the flare are Indicated by a two-character code 
called the importance. The first character, a number from 1 to 4, Indicates 
the apparent area. For areas of less than 1, an "s" is used to designate 
a subflare. The second character Indicates relative brightness; b for 
bright, n for normal, and f for faint. The most recent general discussion 
of solar flares is found In Smith and Smith (1963). 

Hydrogen-a solar flare data are collected by numerous solar patrol 
stations over the globe. The February 1970 Solar-Geophysical Data 
Descriptive Text lists 53 stations. Data, from these stations are found 
tbe Solar-Geophysical Data bulletin and in the lAU Quarterly Bulletin 
on Solar Activity . 

The Solar-Geophysical Data bulletin contains two listings of solar 
flare reports -- those published the month following the flare occurrence 
and those published 6 months after occurrence. The former list includes 
all reported flares of in^ortance greater than one, giving the observatory, 
the date, the start time, the maximum and end times, the solar latitude and 
longitude, the McMath plage region, the duration in minutes, the Importance, 
the observation condition and type, the area, and the maximum width of the 
Bb! intensity. A 1-month listing of subflares is also provided indicating 
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da'te^ time, and location, A full discussion of these parameters is found 
in the Solar-Geophysical Data Descriptive Text which appears in February 
of each year. 

The listing published In Solar-Geophysical Data 6 months after occurrence 
consists of those flare reports filtered according to a program used at the 
Observatoire de 92 Meudon> France, to produce the lAU Quarterly 

Bulletin on Solar Activity flare listing. The format of the 6-month list- 
ing is the same as that of the 1-month flare listing described above. The 
grouped flare data appearing in the 6-month listing are also available 
from WDC-A, Boulder, on punched cards and magnetic tape for the period I 955 
to the present. 

As the number and refinement of observations increases, it becomes 
impressively clear that there is no single, simple flare model but rather a 
large spectrum of flare events which we can only crudely characterize by 
oiir conventional categories of class 1, 2, 3, faint, normal, and bright, or 
by proton, hard x-ray, soft x-ray, white-light, and uv-continuum emissions, 
or by "fast burst", "gradual rise," "fluctuating," and other phenomenological 
features. Many long-held and presumably firmly established concepts of 
ionospheric behavior have become open to question as a result of improved 
observational powers both in space and on the ground. 

It used to be timt ionospheric disturbance phenomena were studied as 
a means of estimating solar-flare x-ray and uv fluxes. But flares have now 
become diagnostic probes for Ionospheric processes because the flare spectrum 
is directly measurable from space and the aeronomic cross sections for ion 
production processes are becoming known with considerable accuracy. 
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Dally Flare Index 

The significant Influence of solar flares on terrestrial disturbances 
of various kinds suggests that a daily flare index, similar to daily sunspot 
number or to the gecmiagnetic activity index, wouH be widely useful. Some 
studies may require construction of a special daily measure, tailored to 
fit the needs of the particular analysis, but a general expression of flare 
activity may suffice for a number of other studies. 

j 

The index published currently at the High Altitude Observatory corres- 
ponds closely to tlie index to be described here. . We shall call the 
present index "the new H.A.O. index," as most of our discussion pertains 
to the original definition. For an individual flare, the integrated Inten- 
sity in Hai is defined as the product of four factors! relative brightness, 
apparent area, a factor depending on the fraction of area at maximum bright- 
ness, and duration. The sum of the Integrated intensities of all the flares 
occurring on a given day is a quantitative estimate of the total excess 
solar energy due to flares. The clear physical interpretation of this in- 
dex is, augmented by other properties that are advantageous in many applica- 
tions. Based on the apparent area and observed duration, this index gives 
a high weight to central flares. Its dependence on flare position is shared 
qualitatively by flare-associated radio and particle emission. The index 
also assigns a high weight to single large flares, as compared to numerous 
minor flares. This weighting also agrees with the probability of occurrence 
and the magnitude of sudden ionospheric disturbances. Thus, this index 
measures a definite physical quantity that is expected to bear a close 
relation to several interesting phenomena. 
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The daily flare index, calculated from the confirmed flares, is 


defined as 






where individual flare area Is measured in square degrees and T* 

Is the effective observing time in minutes, Oni^ those confirmed flares 
of greater than 1 square degree in area, as included in the lAtJ Quarterly 
Bulletin on Solar Activity, are used in calculating the flare index, 
corresponds closely to the flare index developed at the High Altitude 
Observatory to measure the integrated intensity of flare radiation. The 
flare areas sire not corrected for geometric foreshortening, so the 
definition of places great wei^t on large flares, located near the 
center of the sun' s disk. Characteristics of the index are discussed 
in aore detail in the paper by C. Sawyer "A Daily Index of Solar Flare 

p 

Activity" ( J, Geophys. Res . 72 , 385, 1967 ). 

•r 

The table lists the date, index and actual hours of observation 
included in the calculation and follows the table of confirmed Solar 
Flares. 
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Magnetic Condition 

The magnetosphere is disturbed with varying frequency and intensity 
by perturbation in the solar wind. Changes in the solar wind dynamic 
pressure, whether sudden (as in discontinuities) or gradual, result 
in a compression or expansion of the magnetosphere. The effects are trans- 
mitted inward as hydromagnetic waves and are manifested on the ground as 
an increase (for compression) or decrease (for expansion) in the magnetic 
field, representing a magnetospheric response of the simplest kind. 

(l) Magnetospheric Substorms 

Although the exact nature of the interaction between the solar wind 
and the magnetosphere at its boundary is not yet clear, there is increasing 
evidence that particles and other forms of energy are injected into the 
magnetosphere continuously' or in bulk. It appears that when energy accumu- 
lated in the magnetosphere reaches a certain level, it is released rather 
abruptly in what is termed a magnetospheric substorm. The substorm involves 
a large-scale reorganization of the high-energy particle population, of 
the thermal plasma, and of the magnetic and electric fields in the 

i. 

magnetosphere; it also produces a severe magnetic disturbance on the 
surface of the earth at high latitudes, which is called a polar stat or 
magnetic substorm — a complex phenomenon accompanied by other manifesta- 
tions such as particle precipitation, auroral activity, infrasonic wave 
development, and generation of electromagnetic signals. Magnetospheric 
substorms often recur at intervals of 3 or 4 hours, suggesting that this dura- 
tion is the time constant for the buildup and release of energy in the 
magnetosphere. , 
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As Just mentioned, a magnetospheric substorm has a variety of 
manifestations over the entire earth, each of which may provide an 
important clue to the basic processes involved. Recent studies indicate 
that the magnetospheric substorm Is associated with a sadden growth of 
electric fields. These electric fields generate the auroral electrojet 
in the polar ionosphere and accelerate auroral particles, which, in tarn, 
can produce visible auroras and x-ray bursts when they interact with the 
atmosphere. 

The study of polar magnetic substorms alone is not sufficient to 
understand the whole phenomenon of magnetospheric substorms. Intensive 
interdisciplinary studies of polar upper atmospheric disturbance phenomena 
are also needed. Compact, transportable, geophysical observatories 
equipped with identical standard instruments are suggested for such 
studies. A major effort in the analysis of polar upper-atmospheric 
disturbance phenomena should be directed toward the establishment of a 
model, for the distribution of electric fields consistent with various 
mauiifestations of the magnetospheric substorm. The next step is to infer 
how these electric fields are produced as a consequence of , the interactions 
between the “disturbed" solar wind and the magnetosphere. Such studies 
may provide clues to Important problems, including: How can solar wind 
energy enter the magnetosphere? How is this energy stored and finally 
converted into the energy for a magnetospheric substorm? 



(2) Geomagnetic Storms 

The magnetosphere is severely disturbed during magnetic storms, 
which often occur about two days after an intense solar flare, or sometimes 
following the passage of a boundary in the interplanetary plasma sector 
structure. Several hours after the onset of a magnetic storm, a rapid 
succession of severe magnetospheric substorms occurs and the geomagnetic 
field is globally .depressed owing to a formation of a ring current in the 
magnetosphere. Hiysical processes that take place during a magnetic storm 
are extremely complex and many are not well understood. The major pro- 
cesses involved and their effects as observed in the magnetic field at 
the ground are shown schematically in figure 2 (National Academy of Sciences 

1969) . 

A. typical geomagnetic storm has three phases: initial, main, and 

recovery (see Figure 3). A compression of the magnetosphere > manifested 
by a sudden commencement, is the major feature of the initial phase. 

The main phase is characterised by the growth of an Intense ring current 
belt and frequently by severe magnetospheric substorms. During the 
development of the main phase, as many as ten magnetospheric substorms 
may occur in close succession, as illustrated in Figure 3 (National Aca- 
demy of Sciences, I969). These substorras are observed from the ground 
as a repetition of rapid growth followed by decay of the polar magnetic 
substorms. The "AE" index has been, derived to express the activity of 
these polar magnetic substorms and thus of magnetospheric substorms. 



The ring cuirent "belt which develops during the main phase Is pre- 
dominantly coinposed of 5- to 50-kev protons. On the ground, the net 
effect of the ring-current system is a reduction of the horizontal com- 
ponent of the magnetic field In midlatitude and low latitude. An important 
aspect of the ring current belt is a large axial asymmetry during buildup. 
The ”D 4 .” component represents the average intensity of the ring current 
field on the earth* s surface, and "DS" gives the amplitude of the asyme- 
tric. component. There is some indication that the time variation of DS 
is similar to that of the AE index, suggesting that the asymmetric component 
of the ring current is closely associated with magnetospheric substorras. 

,'ki east-west chain of geomagnetic stations is vital for examining the in- 
stantaneous distribution of the asymmetric part of the ring current. 

Typically, about 6 hours after the onset of the main phase, the polar 
magnetic substorm activity begins to decay, while the main phase develops 
for an additional few hours. After the height of the main phase, the field 
begins to recover, sharply for a few hours and then more gradually with 
time. This can taJce several days or even longer. 

In this report the indices D and S are used for magnetically- 
disturbed days and principal magnetic storms, respectively. The magne- 
tically disturbed days (d) are selected in "Solar-Geophysical Data" in 
accordance with the general outline in Terrestrial Magnetism (48, 219, 1943) 
The method in current use calls for ranking the days of a month by their, 
geomagnetic activity as determined from the following three criteria 
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with equal weight: (l) the sum of the eight s; (2) the sum of the 

squares of the eight s; and (3) the greatest K^. (See the discussion of 
the Magnetic Index Ap on page 1-kO for a definition of Kp). 

A table in "Solar-Geophysical Data" presents the principal magnetic 
storms (S) for the month as reported by several observatories through 
cooperation with the International Association of Geomagnetism and Aeronoiny. 



Sudden Ionospheric Disturbances (SID) 

Flares produce a variety of sudden Ionospheric dlsturhances (SID), 
which can he resolved in height from the base of the D region to lower por 
tions of the F region. Timing differences on a scale of a minute or less 
in the. initiation of disturbances and in the rise to peak of the event can 
be highly diagnostic of the natxxre of the ionospheric production and loss 
processes when comparable data for ionizing fluxes are available directly 
from satellites. 

Processes specific to the 60- to 75-km height range are: 

1. Sudden phase anomalies (SPA), where the sky wave changes phase 
with respect to the ground wave as a result of an effective lowering of 
the reflection celling near the base of the D region, sometimes by as much 
as 16 km. 

2 . Sudden enhancement of atmospherics (SEA),, manifested in improved 
reflection of very long radio waves (about 10,000 m) from the bottom of 
the D region. Distant tropical thunderstorms provide a steady background 
on a frequency of approximately 25 kHz, and signal strength may increase 
100 percent during a large flare. 

3. Sudden field -strength anomalies (SFA), observed as interference 
effects over medium distances between sky wave and ground wave when both 
are of nearly equal intensity. As the reflecting celling drops, the two 
waves vary in and out of phase, thus giving large variations in field 
strength at receiver. 

In the 75- to 90-km range there occur sudden ionization Increases 
which lead to: 
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4. Shortwave (5 to 20 MHz) radio fadeout (SWF), which is attributable 
to D region absorption. The phenomena is prompt to within a minute of flare 
outbreak as observed optically froin the earth. 

5* Increased absorption of cosmic radio noise from outer space 
observed by riometers at approximately 19 MHz (SCM) . 

To explore the ionosphere above 90 km we can observe sudden frequency 
deviations (SFD). An SFD is an abrupt increase In the frequency of a high 
frequency radio wave reflected from the F region followed by a slower decay 
to the transmitted frequency. The frequency deviation characteristically 
may exhibit several peaks near maximum and may even show negative deviations 
during the decay phase. In contrast with the D region effects that accom- 
pany 1- to 10-A X-rays, the SFD* s are highly impulsive — the rise to maxi- 
mum is roughly a minute, and the source radiation Is probably in the EOV 
region from 10 to 1030A, most likely dominated by enhancement of the Lyman 
continuum radiation of hydrogen. 

Observations of SFD in radio transmissions via the E and lower P re- 
gions have shown surprisingly fine temporal structure resolvable on a scale 
of only a few seconds. Without knowledge of the solar radiation flux, SFD* s 
could serve as fast-resolving, broadband detectors of explosive flare out- 
bursts, When combined with satellite flare measurements the SFD data can 
provide unique evidence of electron production rates and recombination 
coefficients. Recently it has become possible to record x-ray bursts with 
fine structure similar to SPD’s, and it is remarkable that the SFD and 
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directly measured flux patterns are comparably sharp and can be matched 
to within 1 or 2 sec. From such comparisons it is possible to derive the 
concurrent solar euv continuum emission and to infer ionospheric reaction 
times far shorter than previously suspected. 

SIDs are observed on the sunlit hemisphere and occur almost simul- 
taneously with visual flare observations. SIDs and subsequent recovery of 
the ionosphere have a time duration somewhat longer than flare duration; 
generally from minutes to an hour, with rise more rapid than decay. 

Several of the various SID observations axe routinely identified and 
digitized. , Saraplbs of s^eral types of original records as they appear 
before digitizing are nicely illustrated by Lincoln (I 968 ). The types 
of SID are alphabetically listed below with brief descriptions of each. 

The SCHA, SEA, SES, SFD and SWF listings are available from the WDC-A for 
Upper Atmosphere Geophysics and have also been published in issues of Solar- 
Geophysical Data since 1963. Bursts and SFE are also listed in Solar -G.eophy s 1 cal 
Data . 

Burst - A solar noise burst at rioraeter frequencies is seen as a sudden 
increase of signal strength on the riometer records. This is caused by 
solar radio noise emissions at riometer frequencies. This noise is affected 
by and can mask the ionospheric absorption that is present. The bursts are. 
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of such short duration, however, that absorption changes during the burst 
can be estisnated, in addition to riometer observations, solar noise is 
routinely observed at much higher frequencies than those used for riometer 
receivers. 

Crochet - see SE^. 

SCKA - Sudden Cosmic Noise Absorption, also referred to as CNA 
(Cosmic Noise Absorption) Type I (see below), is characterized on riometer 
recordings by a sudden daytime absorption increase of several db within a 
few minutes. Decay is somewhat slower than rise, and the riometer record 
is perturbed for a period of minutes to over an hour, SCWA involves absorp- 
tion outside of auroral regions. 

SEA - Sudden Enhancement of Atmospherics is observed asan Increase 
in signal strength on wideband equipment operated to detect electromagnetic 
emissions from lightning. This equipment operates in the VLF (10 to 50 kHz), 
range with the most commonly used frequency near 27 kHz. 

SES - Sudden Enhancement of Signal is observed on VLF frequencies 
(15 to 50 hHz). These observations are nearly identical to SEA except that 
the receivers are narrow-band receivers designed to pick up mn-made VLF 
transmissions. As with SEA, signal strength increase is the SID indicator. 

SFD - Sudden Frequency Deviation, caused by ionospheric disturbances, 
is observed by beating the received radio signal with a reference frequency 
that is near the received frequency. This technique allows small frequency 
changes to he observed. SFD is observed in the high-frequency range, 
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and the reference signal is obtained locally from a stable oscillator. 
Occurrence times and deviations in frequency are listed in data presen- 
tations. 


SFE - Solar Flare Effect, or crochet, is observed as a small hook 
on magnetometer records. It is caused by the magnetic field response 
to increased current flow in the E region due to electron enhancement 
induced by flare x-rays, 

SPA - Sudden Phase Anomaly is observed in the same manner as SFD, 
except that VLP frequencies are used and the reference signal is the ground 
wave. For observing SPA on transmitters more distant than possible for 
ground wave propagation (about 200 km), the ground wave reference signal 
Is provided by telephone or a local reference signal source. Is used as in 
SFD observations. 

SWF - Short Wave Fadeout is observed from signal strength records 
of any shortwave (3 to 30 MHz) receiver. Signals from sweep- frequency 
lonosondes (vertical or oblique incidence) may be completely absorbed. 

SWF is categorized as gradual (G), slow (Sl), or sadden (S), depending 
primarily on how rapidly the signal loss occurs. 

CMA 

In addition to SCHA discussed under SID, there are two other types 
of CNA. Data on these latter types are not generally available but are 
mentioned here to broaden perspective. Cosmic Noise Absorption is observed 
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on rlometers that operate between 15 and 6o MHz, but most often near 30 MHz. 
Absorption is normalized to low nighttime absorption occurrences for that 
station. CKA can be divided into BGA (sudden commencement absorption), 

CKA Type I, also called SGNA (sudden cosmic noise absorption), and CNA 
Type II. CM Type II is local nighttime absorption associated with active 
aurora. The absorption may persist for hours with irregular variation 
including individual peaks as high a.s 8 to 10 db persisting for several 
minutes. This isj sometimes called auroral absorption, 

SGA (Sudden Commencement Absorption) is characterized by sudden 
onset, occurrence in zones centered on the auroral zone ( 58 ^ to 75° geo- 
magnetic), and persistence for less than 1 hour. The SCA is similar to 
SCNA in rise time and duration, Type I and II CM are somewhat local and 
irregularly distributed in longitude In contrast to the zonal occurrence 
of SCA. SCAs have been associated with geomagnetic storm sudden commence- 
ment x-ray fluxes. 



2800 MHz (10.7 cm) Solar Flux 

. Within the past 20 years, radio observations, have provided a. new 
means of probing the temperature and electron distributions in the solar 
atmosphere with the resultant possibility of providing other Indices of 
solar aitlvity. In particular. It has been found that the intensity of 
radio emission from the whole disk above a certain minimum value varies 
from day to day and in the vicinity of 3000 MHz the variations correlate 
well with the Zurich relative sunspot lumbers or the total arra of sun- 
spots present on the disk. Sudden enhancements in the emission or bursts 
occur and can be separated readily from the background dally level to pro- 
vide a second radio indication of solar activity closely associated with 
flare appearances. It has been found that the intensity of the daily level 
and of the bursts changes slowly with frequency in the band from 1000 to 
9C00 MHz BO that any observations made on a single frequency in the center 
of this band can be regarded as typical and provides most of the informa- 
tion. The routine observations of the sun at a wavelength , of 10.7 cm (2800 
MHz) were commenced in 194-7 at Ottawa by the Radio and Electrical Engineer- 
ing Division of the National Research Council and have provided observa- 
tions suitable for establishing a radio index. It has been found that 
the single value obtained at local noon is sufficient for most applications. 

The correlation between the orbital deceleration of an artificial 
satellite and the solar radio flux in the decimeter range has led to rela- 
tionships between the atmospheric density at various heights and the 10. 7- 
and 20-cm radio fluxes. Considering that the density variation is due to 
a temperature variation, it can be concluded that the solar extreme ultraviplet 
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radiation (EUV) is responsible for the heating of the upper atmosphere 
and tl^t heat conduction is responsible for its cooling. A change in the 
ultraviolet radiation can be correlated with a change in the radio flux . Just 
as with the sunspot number, which is a well-known index of solar activity. 
The ^0.7“Cm and EUV should be correlated in some degree, inasmuch as many 
features of the sun vary with the solar cycle, including the number of 
sunspots, faculae, filaments, prominences, flares, and height of the chromo- 
sphere and shape of the corona. With regard to the radio waves represented 
and optical manifestations in the EUV regions, how-ever, no exact 
causal relation between these groups of phenomena has been advanced. 

Although it is true that both radiations issue from the upper chromosphere 
and lower corona, our knowledge of the details of their origin is very de- 
ficient. At present, there are no theoretical grounds for an exact relation 
(Anderson, I965). , . ' 

Daily observations of the 2800 MHz radio emissions which originate 
from the solar disk and from any active regions are made at the Algonquin 
Radio Observatory (ARO) of the National Research Council of Canada with a 

i 

reflector 1.8 meters diameter. These are a continuation of observations 
vrfiich commenced in Ottawa in 19^7* Numerical values of flux in the tables 
refer to a single calibration made near local noon at I700 UT. When the flux 
changes rapidly, or there is a burst in progress at that time, the reported 
value is the best estimate of the undisturbed level and provides the 
reference level for measuring the burst intensity. The various types of 
outstanding events are listed separately in another table. The observed 
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flux values have variations resulting from the eccentric orbit of the earth 
In its annual path around the sun. Although these radio values are suitable 
to use with observed ionospheric and other data, an adjustment must be 
Introduced when the observations are used in studies of the absolute or 
intrinsic variation of the solar radio flux. Thus, the tables show both 
the observed flux, S, and the flux adjusted to 1 A.U. , S . The observa- 

Sr 

tions are made for a single North-South polarization but reduced for the 
assumption of two equal orthogonal polarizations. A graph showing the 
monthly highs and lows for the last two sunspot cycles is shown in Fig. 4. 
Relative errors over long periods of time are believed to be + 2 percent, 
over a few days may be + 0,5 percent. The characteristics of the observa- 
tions are surveyed in "Solar Radio Emission at 10. 7 cm" by A.E. Covington 
[J. Ro yal Asl^on, Soc ,, Canada , 63 , 125, 1969 ]» Experiments conducted dur- 
ing the past few years indicate that a multiplying factor of O.906 should 
be applied to the reported flux values in order to derive the absolute flux 


value. 



Heir Solar Flax 

Many workers have carried out the detailed calculations to deterrainf 
how much of the incident solar UV energy is converted to heat in the thermo- 
sphere near 140 km, and how much dissociates and excites oxygen and nitrogen 
moleculf'B. The mean free path in the thermosphere is so large that only a 
small amount of recombination and de-excitation occurs there. In order 
for these processes to occur, the atoms and excited molecules must diffuse 
down into the mesosphere. It is there that this energy is ultimately 
recovered and converted to heat. According to the latest estimates, 
about 70 percent of the solar UV energy is converted to chemical energy. 

Most of the chemical energy is converted to heat in the altitude 
range 70 to 90 km, through tvo recombination processes: One process in- : 

volves the production of Og and 0^* ®ie other process Involves 
de-excitation; 

The extreme ultraviolet line of Hell at 304A contains much of the 
energy involved in heating the atmosphere near 140 km. The measurements 
at 304A were obtained from Timothy and Timothy ‘( 1970 ) for April-August 
1969 , These data are a more accurate index than the 2800 MHz solar flux 
of. the upper atmospheric EUV heating (Anderson, 1965 )* 



Solar Protons (Greater than io Mev) 

Significant numbers of charged particles, mainly electrons and 
protons, are accelerated at the si;n during many solar flares. Some of 
these particles are subsequently detected within the interplanetary medium 
and at the earth. These particles serve as useful indicators of solar, 
terrestrial, and interplanetary conditions. Often, effects from these 
various regions are not separable, but in certain events they are. 

A series of Solar Proton Monitoring Experiments (SIME), initiated 
by scientists at Johns Hopkins University /Applied Physics Laboratory 
and at NASA Goddard Space Flight Center, is intended to provide reliable 
flux and spectral measurements of solar protons over at least half a 
solar cycle and over a wide flux range. Each SIWE consists of an array 
of solid-state detectors that measure the combined counting rate due to 
fluxes of protons and alpha particles in each of the three energy ranges: 
sio Mev/nucleon, ^ 30 Mev/nucleon, and S 60 Mev/nuclepn. The alpha particle 
component is usually very small and may be thought of as introducing an 
ambiguity of less than 10 percent into the proton count rate. Further dis- 
cusslons of this point and of the SI^ in general can be found in the 
"Descriptive Text" of Solar-Geophysical Data (February 1970). 

The SIME was successfully launched aboard the NASA satellites Explorers 
(IMP-F, I 967 - 5 IA) and 41 (IMP-G, 1969-53A) on May 24, I 967 and June 21, 
1969 , respectively. Explorer 34 re-entered the earth's atmosphere on May 3, 
1969 and Explorer 41 re-entered on December 23, 197^‘ Continuous data 
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from the SIME have been presented for the lifetimes of Explorers 34 and 41. 
Similar equipment was launched aboard Explorer 43 (IMP-I, 1971“19A) on 
March l3, 1971 and these data have been used to supplement the Explorer 41 
data, and will continue in 1973* 

The primary purpose of the SPME is to provide systematic monitoring 
of solar cosmic rays over at least half a solar cycle. The basic require- 
ments set forth for such a monitoring program were: 

1) to furnish simple reliable flux and spectral measurements, 

2) to operate over a wide flux range and in particular provide 
coverage for very large events, 

3) to provide a simple and easily reproducible detector system to 
form the basis of an Operational monitoring program. 

All data Cron the SH4E .published as monitoring data satisfy the follow- 
ing requirements: 

1) final orbit determinations have been made, 

2) final time corrections have been'^made, 

3) all data quality flags indicate that the individual data 
point is “good." 

Data that satisfy these requirements are then used to construct hourly 
averages of the respective energy channels. The ^60 Mev and s30 Mev 
channels are sampled for 19.2 seconds once every 2 minutes 43.8 seconds. The 
slO Mev channel is sampled twice every 2 minutes 43.8 seconds at 19.2 
seconds per sample. Thus, the maximum number of points in an hourly average 



is 22 for the a 6o Mev and a 30 Mev channels and 44 for the ^ 10 Mev 
channel. If the number of “good" data points Is less than 5, the hourly 
average is not constructed. 

These hourly averages are then tabulated In a grid where the hours 

of the day run to the right and the days of the month run from top to bottom. 

The day is identified as both day of the month and day of the year. The 

tabular entries are nominally four characters each. The decimal point 

will be, automatically shifted to, the right as the intensity Increases until 

-2 -1 

it finally disappears when, the intensity Is alOOO protons cm sec , If 

- 2-1 - - - . 

the intensity becomes s 10,000 protons cm" sec , the blank apace between 

- 2-1 

columns is utilized so that intensities up ,to 99>999 protons cm sec 
may be accommodated. The entries all beglk with a low intensity format of 

O.XX. The following is an example of how the entries appear fpr an intensity 

-2 -1 ' " -2 -1 
Increasing from 0.55 protons cm, sec to 55>000 protons cm sec . 

0.55 r 

. . 5-50 , f,' ■ . ' 

55.0 

550. 

. 5500 :■ 

55000 

The hourly ayer€ige6 are also plotted vs. time in monthly blocks. 

The response of the ^ 10 Mev channel in the radiation belts is very noticeable 
since this channel also responds to electrons of energy ^ 515 kevl This low 
an electron energy threshold makes this channel respond to the radiation belt 



for a much larger radial extent (at times out to 10 earth radii) than 
chsmnels 1 and 2 whose electron thresholds are several Mev in energy. 

The actual time profiles, due to the radiation "belts, displayed hy the 
^10 Mev channel on the hourly plots are hi^ly variable due to the great 
variations in these trajectories through the trapping regions. 



The Magnetic Index 

lypes of solar-terrestrial correlations and more recent space ob- 
servations demonstrate, when viewed in total, that magnetic activity on 
the earth is very closely related to solar phenomena. This has already 
been discussed to some extent in the section "Magnetic Condition." 

To facilitate correlations with other phenomena, both terrestrial 
and solar, and to have simple criteria for relating changes in magnetic 
activity from hour to hour, day to day, month to month, etc., various 
indices of magnetic activity have been proposed and used. Two of the 
oldest are: C, the daily character figure on a scale of 0 to 2, which is 
a simple numerical, but qualitative, estimate of the day* s activity at 
each observatory, and U which is a measure of the change in the average 
value of the horizon-bal component from day to day. 

In the past two decades considerably more use has been made of an 
index called K, vhen given for a single observatory, and K^, idien given 
for a select group of observatories which we will describe below K 

. pV ■ 

in particular, has been used extensively in space studies both for charac- 
terizing solar-terrestrial conditions and for correlating in detail magnetic 
activity with numerous types of space measurements ranging from magnetospheric 
neutral particle densities to cosmic rays. Thus, it is important that we 
have some understanding of the meaning of K^, 

The K- index of a station Is a quasi-logarithmic measure of the local 
magnetic activity obtained in 3-ho\ir intervals of the Greenwich day. It 
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Is related to the maximum range in the S-hsas- interval of the most disturbed 
component of the magnetic field. The index assumes values between 0 and 
9 corresponding to the quietest and most disturbed conditions, respectively. 
Since the degree of magnetic activity as measured by the range in gammas 
is strongly latitude-dependent, the different magnetic observatories 
must apply appropriate scaling factors to convert the gamma-range to K- 
index in order that the frequency distribution of K-indices be the same 
for all obse37ving stations. Thus, a K of 9 would represent a range of at 
least 300 gammas at low latitudes and 2500 gammas in the auroral zone. 

More than 80 stations currently provide K-indices, 

The planetary index is the average of 3-hourly K-indices from 
12 selected observatories located in middle geomagnetic latitudes from 
47° to 63°. All but one of the stations are located in the northern hemi- 
sphere. The index is coiranonly used as a measure of the global average 
magnetic activity. However, it must be borne in mind that conditions in 
polar equatorial regions may not conform with the planetary K^. 

A linear scale of dist\irbance values a^ has been adopted correspond- 
ing to the quasi-logarithmic values K^. Each Greenwich day has 8 a^-values 
corresponding to the 8 K^-values. The dally average of a^ is called the 
daily planetary amplitude and is denoted A^. High values of this index 
are definitely associated with upper atmosphere heating although the heat- 
ing mechanism is not completely understood (Anderson, 1973). 

The 3-hourly K^ and dally values for each month are published 
(with a 2-month lag) in the Solar-Geophysical Data (Prompt Reports) , 

Part I, CRFL-fPB 135"**; U.S. Dept, of Commerce, Boulder, Colo-, USA 8o302. 
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The AE Index 

Although the Index is capable of describing the. general state of 
planetary geomagnetic activity, it contains contributions from at least 
two major sources, the auroral electro jet and the ring current. To study 
auroral*zone activity it is desirable to maximize the auroral electrojet 
contribution. This has been accomplished by the development of the AE 
index. ®o obtain this index of geomagnetic activity, only slightly sub- 
auroral-zone stations are employed for the most part, and they are chosen 
so as to provide uniformly spaced covereige around the auroral zone. In 
some cases this necessitates the use of southern hemisphere stations to 
take advantage of the apparent conjugacy of magnetic variations of a sub- 
storm nature. 

The AE Index Is constructed by using only the H, component of the 
perturbation field. The H component at each observatory is sealed at 
2 , 5 -min intervals; the average gulet-time baseline is used as a reference 
level. These baseline levels are assumed to be within 10 gamma of the 
undisturbed H component, and thus the actual scaled values of AH include 
contributions in addition to those from the auroral current system. 

All the scaled values from the various observatories are superim- 
posed in a magnetogram format. This composite drawing is enclosed by upper 
and lower envelopes representing the maximum positive and negative values 
of AH for all observatories at each given time. The amplitude of the 
upper envelope at any instant is denoted by AU, and the amplitude of the 
lower envelope is designated by AL. If the station distribution were 
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closely enough knit, AU and AL would represent the maximum positive and ne 
gative deviations occurring along the auroral' zone. 

Physically speaking, AU gives a good representation of the maximum 
niagnetlc perturbation generated by the eastward electrojet usually found 
In the afternoon sector. Similarly, AL represents the maximum magnetic 
perturbation generated by the westward electrojet in the morning and 
midnight sectors. At the present time it appears as though the eastward 
and westward electrojets may fluctuate Independently of one another, 
and it may be useful to treat AU and AL as independent indices. However, 
at the preseht time it is customary to combine these indices to give a 
direct measure of the total njaxlmum amplitude of the' eastward and west- 
ward electrojet currents. Hence, the index AE is defined by AE = AU - AL. 
Thus, AE represents the difference in levels (measured in gammas) between, 
the upper and lower envelopes at any given instant in time. 

Although AE is generally calculated at 2.5-min intervals and is 
s^Q’ilu^ls in this form from the World Data Center, hourly average values of 
AE are also available for the period 1957-196^, inclusive (University of 
Alaska Reports UAGR-I 92 , UAGR-200). More recent AE index values through 
1970 are available from the National Space Science Data Center, Goddard 
Space Flight Center, Greenbelt, Md. 20771 . 

Observations of traveling ionospheric disturbances emanating from 
polar regions have led to the suggestion that auroral energy may be coupled 
to atmospheric wave motions through heating produced by ohmic dissipation 
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of electric currents, Joule heating. Cole’s (I962) pioneer work has 
established that Joule heating in the region 100-200 km could lead to 
large temperature chaxiges, through heat flux by thermal conduction in 
the vicinity of I50 km and above. However, below I50 km it appears that 
the relaxation time, i.e., the time it takes an initial temperature impuls 
to fall to e is at least one day and is independent of T^. On 

this basis, short period gravity waves (---- 1 hr) should transfer energy 
more efficiently than conduction below about I50 km, if the vertical 
scale of the propagating waves is approximately greater than the local 
scale height and if the coupling of auroral heating into gravity wave 
energy is efficient. Blumen €ind Hendl (I969) have investigated Joule 
heating in the vicinity of the auroral electro jet (IOO-I50 km). 



The D . index 
— st 

The index provides a planetary magnetic index on a hourly basis 
that has been used extensively in studies of geophysical disturbances, 
particles and plasma in the magnetosphere, solar^terrestrlal relationships, 
and cosmic rays (Suguira and Cain, I969). 

Equatorial is a measure of the mean departure from normal of 
the horizontal component, H, of the Earth* s magnetic field observed at a 
group of low-latitude stations, whereas is based on 3 hr ranges of the 
field at stations in higher latitudes. 

The strength of lies in Its ability to detect all magnetic 

storms. This stems from the fact that the ring current is world wide in 

nature, and thus there is little or no possibility, of D . not reacting to 

the growth of the ring current. Furthermore, within about two hours, the 

index should be able to identify the onset and termination of the main 

phase of a magnetic storm (during which time the major portion of the energy 

of the storm is dissipated in the magnetosphere). Thus, the D index 

sx» 

gives a good qualitative description of the gross level of magnetospheric 
activity at any time, although it cannot be used to reveal the presence of 
individual substorms. 

The values of used in this analysis came from Goddard Space 
Flight Center reports (see for example, Sugiura and Poros, I971). 



Solar X-Rays (1-8A) 

This index is used to indicate the magnitude of the x-ray flux 
between 1-8A -irtiich can cause atmospheric heating between 55-105 km. Solar 
x-ray emissions, like solar radio emissions, originate from both thermal 
and nonthennal processes that take place primarily in the solar corona. 

Solar x-ray emissions can be divided into three components: the quiet 

sun component, the slowly varying component, and the burst component. 
Nonthermal x-ray emissions occur only in the hard x-ray emissions (photon 
energies greater than about 10 kev or wavelengths less than about lA) 
accompanying solar flares. Solar x-ray astronomy has been reviewed by Krimi 
gis and Wende (I 97 O). 

Solar x-rays are directly observed with satellite -born detectors 
and indirectly observed through their effects on the ionosphere. The 
occurrences of solar x-ray bursts are Indicated by the sudden ionospheric 
distui-bances (SID). Patrol types of solar x-ray measurements must be 
carried out using instruments flown on satellites. Because solar x-ray 
instrumentation must be flown on satellites, the data quite often contain 
periodic and signifiant time intervails when no data were obtained due to 
the satellite’ s being eclipsed. 

Ion chambers, proportional counters, Geiger tubes, and^ scintillators 
are used to measure the integral solar x-ray flux within specific wavelength 
or energy passbands. A review of such instrumentation is given by Weupert 

(1969). 
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Solar x-ray fluxes, integrated over the entire solar disk in units 

-2 -1 

of ergs cm sec over a specified passband, have been obtained from 
several satellites. The x-ray fluxes are useful and sensitive indicators 
of the level of solar activity. 

The data used in this report consist of tables of hourly averages 
of x-ray flux in the 1 to BA band. The primary source of the data is the 
Naval Research Laboratory's SOLEiAD 9 and 10 satellites. Each satellite 
stores 0.5 bo 3A, 1 to 8 A, and 8 to 20A data in the satellite memory with 
a one minute time resolution. Therefore, a continuous record of the x-ray 
emission from the sun, except for gaps due to satellite night and 
charged particle interference, is available for these three bands. A 
complete description of the SOLRAD 10 experiments is given in NRL Report 
No. 7408 titled "The SOLRAD 10 Satellite, .Explorer 44, I 97 I-O 58 A." A 
complete description of the SOLRAD 9 experiments is given in KRL Report 
No. 6800 titled "The NHL SOLRAD 9 Satellite, Solar Explorer B, I 968 -I 7 A." 

The data from the 1 to 8A ionization chamber are converted to a 1 
to 8A energy flux based on a 2 x lO^K gray-body solar emission spectrum 
(Kreplin, R.W. , Annales de Geophysique , 17> 151^ I 961 ) . The averages 
include data obtained during solar flares but data contaminated by charged, 
particle interference are excluded wherever possible. When SOLRAD 10 
data for an hour are missing, SOLRAD 9 data are used. Measurements by the 
SOLEIAD 9 and SOniAD 10 1 to 8A ionization chambers are virtually identical 
so normalization is not re<iuired. For identification purposes only, the 
SOLRAD 9 data are followed by the symbol <. The v^ues given for each 
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hourly average are in units of lo'^ erg^/cm^ sec, and the time scale Is 
UT. The-tabular entries are made with a shifting declmal'point. There 
fore, the values which may be expressed in the four spaces assigned to 
each entry range from .001 to 9999 times the basic unit, lo‘^ ergs/cm^ i 
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5. SOLAR-TERRESTRIAL LINKS AND STRATOSPHERIC WARMINGS 

Before physical theories can be invented, to explain certain suspected 
solar-weather relations, including stratospheric warmings, several links 
must be filled in. Some of the principal missing links in the development 
of such theories are as follows: 

Thorough knowledge of incident solar flux - Long-term accurate measure- 
ments are needed at all wavelengths, including Lyman-alpha (1216A), EUV (lOO- 

90oi), and x-rays (< lOOA). Most of the Lyman-alpha is absorbed below 100 km. 

-3 -3 

Its flux varies from about 6 to 3 ergs cm sec between high to low solar 
activity. Changes during flare periods are not significant geophysically. 

The x-rays should be measured, globally, since the greatest variation in 
solar photon output occurs in the x-ray region (up to 3 orders of magnitude 
with periods of minutes to hours). It is necessary, to separate out the 
ionizing effect of the x-rays from that of longer wavelengths. The main 
action of x-rays is to create ionization between 120 km (lOOA) and 55 km 
(lA). However, they are not the sole Ionising agent in this altitude range 
nor is their contribution to heating the major one. 

Appropriate data on the primary photon interactions - Quantitative 
data on the rate coefficients of various types of transitions for all the 
Important constituents at all important wavelengths so that the consequential 
effects of the primary photon interactions, such as ionizations of atoms 
and molecules, dissociation of molecules, etc,, can be assessed. 
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Heat budget of the mesosphere near 80 km - What are the altitude, 
latil.ude, temporal, diurnal, and ceasonal variations of* mesospheric tempera- 
tures? How are the temperatures related to the intensity and variai;ions 
or the solar layman-alpha and x-rays? Are there other important heat sources 
for the mesosphere, such as particle Influx, and how do they contribute to 
the temperature structure? What is the energy balance of the mesosphere as 
a function of height, and what are its global variations? A difficult aspect 
of this problem is the transfer of energy as a result of radiative processes 
and various circulation processes. The altitude at which energy is absorbed 
need not in general correspond to the altitude at which it is converted into 
thermal energy. 

Circulation In the Mesosphere ' - It is necessary that the temperature, 
density or pressure, and composition of the mesosphere be measured as func- , 
tlons of altitude, location, and time, In sufficient detail to disclose the 
general patterns and systematic variations that occur with (a) time, 

(b) season, (c) sun-spot cycle, (d) location, and ( 5 ) the condition of other 
parameters (meteorological, auroral, geomagnetic, etc. ) . The limits and 
characteristics of randcan variations that cannot be associated with other 
geophysical parameters should also be determined. Thus, it is necessary 
to establish the inherent variability of the mesosphere, both short term 
(to determine possible effects due to x-rays) and long terra (to find possible 
effects due to Lyman-alpha). 
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What are the vesrbical distributions of the chemically active atoio - 
spheric species ? - What processes control those distributions? The most 
important species are ozone, nitric oxide, carbon dioxide, hydrogen and its 
compounds, atomic oxygen, argon, methane, and excited atoms and molecules. 
These minor constituents need to be measured as function of latitude and 
time. Water vapor and methane dissociate in the 6O-IOO kra region, and the 
dissociation and recombination rates are so poorly known that a direct 
measurement of their concentration is needed. Accurate measurements of the 
vertical profile of some or all of these species will be interpretable in 
terms of the large-scale circulation in the mesosphere. 

Laboratory measurements of reaction rate coefficients - While some ■ 
of the chemistry of the mesosphere is understood, there are serious gaps in 
our knowledge of. the chemical processes involved, and many of the needed 
rate coefficients are missing. 


How is turbulence generated in the mesosphere ? - If waves are respon- 
sible, what are their origins, and are their amplitudes sufficient to 
cause nonlinearities? 

How does eddy diffusion affect the distribution of atmospheric gases 
and the transport of heat ? - This is an important question near 80 km. It 
is necessary to determine the value of the eddy diffusion coefficient as a 
function of height and time from 6O-IOO km. 
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jtot _are the physical mechanism(s) responsible for interactlon( s) 
between the mesosphere and lower atmosphere (troposphere and Gt“ratosphere)? - 
How does their relative importance vary with aJ-titude, location and tiiuo? 
Statistical studies indicating correlation between upper and lower atmospheri 
phencanena have little validity unless backed up by physical mechanisms. For 
example, the statistical correlation between the international magnetic 
character index C (which is coirelated with solar activity) and lower atmo- 
sphere circulation systems may result from dynamic interconnections between 
lower and upper atmosphere circulation systems rather than of solar emission. 

Hydrodynamic theory of atmospheric perturbations - This theory needs 
to be developed, in order to find the Influence of upper atmospheric solar 
heating, etc., upon the circulation in the atmosphere. How is the general 
circulation influenced by irregular changes in the sun's emlssivlty or by 
variations of the reflectivity of the earth-atmosphere system? The physical 
basis for changes in lower atmosphere circulation patterns is not fuliy ■ 
understood. 

Development of techniques to evaluate the significance of weak 
relationships or relationships exposed from a small number of cas es - 
Only about .01 percent of the solar radiation is in the wavelengths con- 
sidered. Even if this energy absorbed near 80 km does influence the earth' s 
lower atmosphere, how could this be detected? 

Figure 5 indicates the many types of perturbations affecting the 
stratosphere and the connections among them, . 
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52 
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V' - WEST GKIPviNISH (57.\, 7\n ROCKET DATA 


— 1. 7 DEC 67 -BEGINMXG PHASE. MAX WARMING (60 KM) 
IS -fire (CURVE FOLLOWS GROVES CURVE (5) FROM 
.ABOUT 35 KM TO 55 KiM). . . 

2. 13 DEC. 67 -^LVX. WARMING HAS DESCENDED TO 

■46 KM AND INCREASED TO +18“ C. 

3, 21 DEC, 67 " 1 VCREASING PHASE MAX. WARMING 

(43 KM). 15+40“ C 

4. 27 DEC. 67 - DECREASING PHASE MAX, WARMING 

(44 KM) DECREASED TO -4®C. 

5, GROVES (1S70) AVERAGE ATM. 60N DEC. 1 


TEMPERATURE CC) 

FIGURE 1 WEST GEIRINISH ROCKET DATA DURING WARMING (Latitzke 1971) 
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FIGURE 3 An example of a typical geomagnetic storm showing 
the three characteristic phases and the relationship between 
the Dst, AE, and DS indices. The Dst index is a measure of 
the magnetic field intensity of the ring current and also of 
the total kinetic energy of the ring current belt. The AE 
index is a measure of the intensity of magnetospheric sub- 
storms. The DS index is a measure of the magnitude of the 
asymmetry of the ring current field. 
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CURVC <-- --lAPP«OXIMATeLV SEPARATES SUNSPOT COMPONENT FROM SASIC COMPONENT. 
^MACNITUOE OF THE RADIO QUIET SUN AT SUNSPOT MINIMUM SHOWN 

FIGURE 4 Monthly Maxima and Minima of 3800 MH^ (lO-T cm) Solar Flux, 1947-I972. 

(Covington, 1969) 
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FIGURE 5 Perturbations Affecting the Stratosphere 
(Systems Control, Inc,, 1973) ' 
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1. WAFMIHG CASE STUDIES 


The study began vith the selection of a set of cases for study. The 
original set of cases vere identified throijigh an extensive survey of the 
Uterature on Stratospheric Warmings. Through detailed study of these cases 
and additional data acquired during the study, the case histories were expanded. 

A number of data sources were used in the study of each case. Routine 10-mb 
analyses on a dally basis are available In mlcroflLn form for the northern 
hemisphere from the Environmental Data Service, National Weather Records Center. 
These maps were augmented by a number of additional analyses of warmings kindly 
leaned the Investigators by Mr. Roderick S. Quiroz, Upper Air Branch, National 
Meteorological Center. Ebctensive use was made of published materials from 
Journals. These data were derived, from rocket sondes, balloon sondes, and 
satellite analytic soundings and radiance measurementa. The bulk of the data 
used for southern hemisphere analysis was loaned to the Investigators by 
Dr. Sigmund Fritz, National Environmental SatelUte Service and consisted of 
daily digital global maps of SIRS radiance. This data was augmented by a small 
number of Nimbus 3 MRIR plots of radiance from Channel 3 (l^.5-15.5^m) obtained 
from National Space Science Data Center. Attends to obtain vertical tenq)erature 
profiles through an identified wanning from the Environmental Data Service vere 
unsuccessful. Also unsuccessful were attenq>t8 to obtain any data from Clarendon 
Laboratory, University of Oxford, the principal investigators of the Selective 
Chopper Radiometer on Nimbus. 

The Stratospheric Warming of April - May 1969 

Fritz (19Y0) describes a stratospheric warming present in early May over 
the Indian Ocean neax 45“S. Figure 1 presents a march of the maximum zonal radiance 
(sms A, Channel 8) for kO‘S, JO'S, 6o*S, JO’S, and 80”S from 15 April - ik May I969. 

II-l 



This data was derlvted from digital plots of SIRS radiance obtained from 
Dr. Sigmund Fritz, National. Environmental Satellite Service. The warming is 
evident in early May from 1K)-60*S. However, from 70-8o“S, the maximum radiance 
occurs on 2 k April. Figure 2 displays the maximum radiance at latitudes of 
^0®S or greater. Five peaks are apparent; the two major peaks occur at 
30 April and 5 May. Unfortunately, the lack of radiance data for days earlier 
than 13 April does not allow an accurate selection of the date of the start of 
this vanning. Therefore, no detailed analysis will be attempted; however, 
some of the data fran this warming will be used in later general analysis. 

The Stratospheric Warmings of July August 1969 

Three distinct stratospheric warmings occurred between 3 July I969 and 
19 August 1969^ "With peaks near I3 July, 6 August, and I6 August (Pigs, 3 and 4). 
With regard to the July wanning. Pig. 5 (Fritz and Soules, 1970), displays the 
change in radiance in the COg band centered at 669*3 ^ measured by the Nimbus 3 

SIRS instrument from 25 June - 10 July I969. Fig. 6 (Wark, I970) indicates that 
this band (channel 8) is most representative of temperatures near 10 mb (31 km). 

In Pig. 5, cooling (c) occurred everywhere from 20**3-4o®N, vdaile large-scale 
wanning took place over much of the region south of 20‘*S. The 669 cm"^ radiance 
field on 8 and 9 July I969 is presented in Finger and Mclnturff (1970), A 30-mb 
{ 2 k km) temperature field, for 9 July I969 is shown in Miller et al. (l^O). 

The starting and ending dates of the July warming may be determined from 

Fig. 3, ^ere the soUd curve la derived from the daily maximum radiance (ergs 
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cm s ster cm ) values at latitudes 140®S or greater from SIRS A, Channel 8 
world-wide digital printouts, displaying radiance values at intervals of 10 degrees 
of latitude and 30 degrees of longitude. The maximum radiance values were between 
40-60 *S latitude. The dashed curves represent the 10-mb (3I km), 30-mb (24 km) 
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temperatures from Mirny, Antarctica (Finger and Mclnturff, 1969; Smith et al., 
1970). From Fig. 3 the start of the warming at 3I km is (solid curve) at 5 July; 
the maximum at I3 July, and the end at 21 July. The 10 -mb and 30-mh curves indicate 
that the warming started 3 Jiay. The 10 -mh curve indicates that the peak of the 
warming near the South Pole at 3I km was near I6-I7 July. Close agreement Is not 
expected hetveen the solid and dashed curves since the respective locations differ 
appreciably. 

From Pig, k, solid curve, the second warming starts on 1 August, peaks on 
6 August, and ends on 11 August. The thir^ warming starts on 11 August, peaks 
ends on I9 August. For these latter two varmings, the Mimy 
10-mb and 30-mb temperatures agree within about one day for the start, peaks and 
ends of the vaimlngs. Fig. 7 presents a 30 -mb ten^erature field for 17 August I969, 
near the peak of the third warming (Miller et al., I970). Two warm centers are 


shown. 



Correlation of Geophy sical Indices with SIRS A. Chatmel 8 Radiances , July - 
August 1969 . - 

Table 1 lists the Indices vs, day of the month from I-I 7 July I 969 . 

No daily flare index was available for this period. If the start of the 
vanning was 3 July (lO-mb and 30-mb curves. Fig, 3 ), then possible heating 
candidates are solar x-rays and/or electrojet heating. Solar flare of impor- 
tance 2 or greater occurred on July 1, 3, 4, 5 , 6 , and 7 . On 3 July, there 
vas an SID of Importance 3 and a strong x-ray flux (1-8 a). Table 2 shows 

the 1-8A solar x-ray bursts (Explorer 37 ) from I-I 6 July I 969 . The 3 July 

-2 -1 

flux of 0,19 ergs cm s was not. only the strongest x~ray burst measured 
during this period but It also was of long duration (40 minutes). The AE 
maximum of 294 on 1 July Indicates that joule heating associated -with the auroral 
electrojet must be considered also as a trigger mechanism. On 5 July, there 
vas a flare of Importance 3, SID of importance 2"^, and maxima in solar x-rays, 

He li flux and AE. If the heating started on 5 July, as. indicated by the 
solid curve of Fig. 3, ; then EUV' heating (He II flux) (140 km) and electro jet. 
(AE) (100 km) are possibilities. Fig. 3 (Solid curve) indicates abroad 
warming maximum from 8 -I 5 JiHy. This broad maximum could result from enhance- 
ments indicated by indices for x-rays. He II flux, Ap, D , and AE, The 
x-ray fluxes (Table -2 ) do not appear strong enough during this period. 

Again, electro jet heating offers the best possibility since it occurs lower 
(100 km) than the EUV heating (140 km). 

Table 3 lists the indices vs. day of the month from 29 July - 
14 August 1969 . i'ig. 4, solid curve, indicates that the start of the 
second warming was 1 August; it was' followed by a very steep rise in radiance 



to the peak on 6 August. A solar flare and SID of Importance 3 occurred on 
29 July. Also, an x-ray naxlnam occurred on this date, but the xlray bursts 
were not very strong. On 30 July, a principal nBgnetic storm occurred to- 
gether with maxima in Ap and AE. On 31 July, a Dst naxlmum occurred. The 

most likely cause of the start of the second warming was electrojet (AE) 
heating. 


The third vanning started on 11 August and peaked on l6 August 
(Fig. % ). Both, the solar x-rays and He II flux had maxima on 11 August. 
However, the associated fluxes were not particularly strong. On 7 August, 
there was an importance a solar flare and a weak maximum in x-rays. On 9 
August, Ap had a maximum. On 10 August, Dst and AE had moderate maxima. 



The Stratospheric Vfarming of December I 969 - Jemuary I 970 

A mid -winter warming affecting a deep layer of the Northern Hemi- 
sphere stratosphere occurred in December I 969 - January I 97 O. While limited 
circulation changes were observed, meteorological rocketsonde data indicated 
a warming of significant proportions, with temperature Increases of up to 
80°C near 40 km. At the 10-mb level (approx. 3l km), unusually warm tem- 
peratures were first noted over southern Europe after mid-December, Accord- 
ing to Gelman ( 197 O), the warm area expanded and the center moved norths 
eastward along the direction of wind flow to central Siberia, where temperatures 
warmer than - 10 °C were observed, at 10 mb around S7 , December, ..At West Geirinish, 
Scotland, extreme conditions were observed, with a temperature , of +38°C and 
winds in excess of 300 kt, at 43 km« • An Intense anticyclone centered over 
Alaska had developed during the course of the warming. At the end of Decem- 
ber the warm air at 10 mb was displaced poleward from center Siberia in 
the southerly flow between this anticyclone and the elongated polar vortex., 

By January 3, the warm air had readied the polar region with the moderated 
cold air being displaced southward, thus reversing the normal wintertime 
temperature gradient. 

Miller et al., (1972) constructed high-altitude analyses of the warm- 
ing event, up to the 2-mb level (43 km). The Satellite Infrared Spectrometer 
(SIES) aboard the Nimhus 3 satellite yielded large amounts of radiation 
temperature data (Quiroz, I 971 ). The use of rocketsonde and satellite data 
combined with high-altitude radiosonde observations allowed hemispheric 
synoptic analyses at the 10-, 5-j arid 2-mb levels to be constructed. 
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Unfortunately, it was not possible to construct the high-altitude analyses 
on a daily basis. Consequently, analyses were done for the seven, selected 
days: December 17 , 24, 29 , 3l and January 3, 7 1^- 

Figure 8 depicts the height and temperature analyses at 2-mb for 
24 December. At this time the polar cyclone was the prominent . circulation 
feature with ridge development near Japan. The cold centre associated 
with the polar low is located over Canada, with warm air centers over the 
USSR and southeast United States. Comparison of the positions of these 
warm. and cpld areats with those at ^-vib (not shown here) indicates a general 
westward slope with height. 

By 29 December, Fig. 9? the Eastern Hemisphere warm air center had 
intensified considerably and moved westward with time. As an example of the 
change in Intensity, a rocket observation at the United Kingdom's West 
Geirinish station on 23 December measured a 2-mb temperature of -23*^C; 
but on the 27 th, the observed temperature was +23^C, a rise of 46^C. In 
addition, the general height pattern developed a more wave 2 type configura- 
tion, with a general elongation of the polar cyclone. 

Figure 10 illustrates the analyses for 3 January. It is obvious that , 
tremendous changes have taken place in the height field in the 5-day period 
from 29 December. The most sbrlklng feature on this chart is the segmentation 
of the height field. Remnants of the. polar cyclone were located mainly 
over the Northern Pacific and Atlantic Ocean areas, i^ile relative highs 



were situated over the North American and Eurasian continents, including a 
ridge across the polar regions. The temperature fields had continued their 
westward retrogression with the warm air, about +20°G, mainly over North 
America and the cold cell associated with the deepest cyclonic segment 
located with the coldest air over the Alaskan-Slberian area. This series 
of maps clearly points out the rapidity with which the wantilng and circu- 
lation breakdown took place. 

The last analysis of this series. Fig. ll, shows that by the middle 
of J^uary, the warming in northern latitudes had run its course and cool- 
ii^ was taking place. The height field, however, was still ' considerably 
perturbed and the wintertime polar cyclone had not yet returned to the Arctic 
region. The circulation changes In lower, levels of the stratosphere (at 
and below about 30 mb) were not nearly as pronounced as those described 
above. Although warn air was noticeable over the Siberian-Alaskan region, 
the polar cyclone remained strong and primarily in its seasonal Arctic loca- 
tion. 

Figure 12 shows the change in temperature at 2 mb (4-3 km) from 24 De- 
cember 1969 to 3 January 1970- This map Was derived by subtracting the 
temperature field in Fig. 8 from that of Fig. iQ. The maximum temperature 
increase is 80 degrees with center located at 258 longitude. The ni£txxmum 
decrease Is around 30° with center at l 82 ° longitude. 



A description of the major radiance changes during December I969 - 
January 1970 has been given by Quiroz (1971 )- The analysis was based on Nimbus 3 
SIRS maps. The features of greatest Interest Include: 

1) Development of a strong bl*cellular pattern of high radiance with 
eastward displacement of the system centered over northwest Africa. 

2) ^parent wave reinforcement over Siberia (Fig. 13) resulting in a 
radiance maximum > 70 ergs/ in channel 8, around 27-28 December. 

3) Translation of high radiance system across the polar area to the western 
hemisphere, with maximum local increase of 38 ergs/ near 75*N, i3e, from 2k December 
to k January. 

4) Daii5)ening of this system, with generation of westward-traveling waves 
across North America to Asia. 

5) Displacement of polar lov-radlance system (central value 31,5 ergs/ 

on 2U December (Fig. I3) to Canada, thence westward to Asia (Fig. l4) and return 
to pole by 9 January. 

The down-pointing NIMBUS SIRS instrument has a field of view of about 200 km 
by 200 km. (The angle of view Is i6*,) The orbit Is sun-synchronous and is 
therefore always shifting westward as the satellite proceeds northward to 8l*N, 
then 'southward on the opposite side of the globe to 8l®S. The longitude separation 
between successive orbits is 27®. The instrumental accuracy is high. Observations 
were nude every 8 sec along the orbit, yielding some ten measurements every 4® of 
latitude. The rebdleuice data which were used are generally averages of up to five 
of these measurements, taken within squares defined in the numerical prediction 
grid used by the National Meteorological Center, about 38O km on the side. All 
orbital data that, met basic calibration criteria, taken in a 24-hr period centered 
on 1200 QfT, were used. 
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printout vas obtained at 2-day IntervalB In a 6o-day period centered 
on 26 December I 969 near the time of maxlimim vnrmlng. Analyzed channel 8 
radiance maps for tvo selected dates are shovn in Figs. 13 and l4. 

Tor a typical varm-layer thlotess (20-25 km), Quiroz (19TI) has shovn 
that the central altitude and the anqplltude of warming can be satisfactorily 
discriminated through the use of two sinqple parameters Including the ratio of 
, radiance change In NIMBUS SIRS cftiann 7 {-678 cm'^) and 8 (-669 cm’"^). 

A nomogram is developed for speclfylaag warming altitude and an5>lltude from 
observed radiance changes and is tested during the December 1969 - January 1970 
wanning event. 

A distinct pattern is revealed In the ratio of change In channels 7 and 8, 
with decreasing values poleward for the period 20-28 December and from the pole 
to the NorUi Atlantic in the next 8-day period. 

The ratio viaries from more than 1.00 in the "rear" portion of the system 
^ large Increase, to less than 0.30 In the forward part. On the basis of the 
, Quiroz model, this ratio pattern together with the actual radiance Increases 
la channel 8 implies warming of approximately 6o"C centered near an altitude 
of 25 km In the rear portion; and of 70-80*C centered Just above 4o km in the 
forward portion. 



. Northern 10-mb (3l tan) chains were available from the 

, Hi^erical Weath^ Prediction Unit, National Weather Service, for 1200 GMT 
ea^^d^ ^ini)ecen^ I 969 ^d Jahiilary l 970 .These series of maps offer the 
■bestrop^rt^ity define ipany of the SW characteristics in detail. Inas^ 

; much as the , war^ proceed from above, these charts can be used to 
determine the latest dates at which the warmings start. Also, they probably 
give the best Information on the end of the warmings. Figure I 5 shows the 
chart for 26 December 1969 ; the maximum warming center (266^K) is located 
at 100°E and -57^* Figure I 6 illustrates, the movement of this center from 
/23 Dec^her :(l4o°E, ., 50 °W)/to 30 -I)eceiriber Several other warm- 

ing centers were also present at 10 mb during the same period, but thfey were 
not as intense. The circles on Fig« I 6 indicate the location of warming 
centers at 2 mb (43 tan) for 24 and 29 December I 969 . It is obvious that 
the maximum warming centers on these dates do' not coincide at 10 mb and 2 mb. 
Figure 17 is a plot of; the daily maxlmuin temperature at 10 mb (3l km) 
occurring in the higher latiti^es of the Northern Hemisphere, The maximum 
temperatures on 21 December I 969 and 7 January 1970 the same, as indi- 
cated by the straight line. It appears that the start of the warming occurred 
on about the 17 December and the end by the 5 of January at 31 tan altitude. 
The dashed line on Fig, I 7 ^indicating the maximum 2 mb (43 km) temperature, 
is b a ;a e d on 5 values from Fig.8-11 and Miller et al. ( 1972 ) (indica-. 
ted 'by crosses. on Fig. 17). It appears that the highest maximum temperature 
at 2 mb occurred on Dec. 27,.^ere was little change between 29 Dec. I 969 
and 3 Jan. 1970> diuring the same period ■vrtieii the lO mb maximum temperatures 
are decreasing markedly. Although some of the difference occurring between 



10 mb and 2 mb may be explained by the scarcity of 2 mb data, it is diffi- 
cult to explain all the difference on this basis. Rather, it appears likeOy 
that warming center locations (Fig. l6 ) and msignltude (Fig. I7 ) differ 
considerably at the two altitudes for the same date. 

Hence, it can be concluded from 2 mb ,’ 10 mb, and radiance data 
(Quiroz, 1971) that the maximum warming occurred between 26-29 December 

1969 at altitudes 25-^3 km. It is more difficult to determine the start 
or . end of the warming. The warming certainly has started by 24 December, 
although the 10 mb data (Fig. I7.) Indicates that it could have started as 
early as I7 December. Likewise, the end is difficult to pinpoint. Figure 
11 . and Fig, ,17 indicate that the wanning had run its course by 14 Jan. 

1970 at 2 mbj it definitely appears to end at 10 mb by 5 Jan. (Fig, 17.). 
Rothiog can be stated concerning the altitudes above km. 



Correlation of Geophysical Indices with 10-MB Temperature Increases , 

December 1969 

I 

Maxima in the dally flare index, 8XD* s and x-rays are expected 
to occur on the same day as solar flares. This is not necessarily true 
of the 2800 MHz flux unless the EUV associated with the flare is very strong, 
Llkevlse, maxima of solar protons, A^, ©st, and AE should occur about two 
days later than a flare. This Information may give clues that can be used 
to distinguish among possible heating from the various solar effects. For 
example, we should expect x-ray heating to occur immediately after a solar 
flare, while proton heating would occur approxlmateOy two days later. Also, 
x-ray heating (electromagnetic) should occur on the dayside of the earth 
only. Protons precipitate on both the day and night sides with a marked 
preference for the auroral zones. 

Table ^ lists the indices vs. day of the month from II-27 December 
1969. On 17 December, the 10-mb (31 km) maximum temperature was It 

increased 10 degrees to 2i6“K on the l8th of December and dropped down to 
2^a'k: on 19 December. The warming center was located near 35*H, 20^< Accord- 
ing to Groves (1970), the average temperature at 35*N during December is 
about 237*lt« Assu m ing for the moment that this marked temperature rise was 
due to an extra terrestrial (solar) cause, ^diat is a possible candidate? 

Note in Table that an importance 3 solar flare occurred on 
14 December. This flare occurred at 0345 UT in McMath plage region 10477. 

On this same day the dally flaj:e Index and solar x-rays (I-8A) ej^rlenced 
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relative maxima of II5 and 9 x 10 ergs cm s , respectively. About 
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two days later (Dec. l6), the solar protons greater than 10 MEV exhibited 

-2 -1 

an hour3^ maximum of 159^ cm s'. Also, relative maxima were exhibited on 
l6 December for Dst, and AE. On 17 December, an importance 2 flare occurred; 
the dally flare index was 457 or almost four times as high as on 14 December. 
Also,an SID of importance 2 and solar x-ray maximum (26 x 10“^ ergs cm"^s~^) 
occurred on I7 December. The proton maxima> recurring at 2-3 day 

intervals, are probably due to the earth^ s radiation belt. 

Therefore, if the flare of December 17 was responsible for 

the temperature maximum of 248K occurring on the 18th of 

December, the most likely mechanism was heating of the atmosphere below 100 km 
by‘ solar x-rays (1-8 a) occurring on I7 December. Table f indicates that 
the x-rays on I7 December were Inore intense and longer lasting than on 14 Decem- 
ber. If this were true, then it took about a day for the heating between 
60-100 km to reach the 10-mb surface at 3l km. 

If the heating is associated with the flare of 14 December, then the 
AE index maximum of I90 on I6 December indicates that the most likely avenue 

was Joule heating associated with the auroral electrojet near 100 km. The maximum 
hourly value of the AE index (52l) occurred on I6 December at I7OO UT. 

The AE index was above the December 1969 monthly mean of 110, 71 percent of 
the hours on December 16 (Table -6 ). Thus, if electro jet heating was the 

cause, it took about two days for the heating to reach the 10-mb surface 
at .31 km. One problem with electro jet heating is that it is restricted to 
higher latitudes near the auroral zone, at least . initially. However, the 
warming center on the I8 December was located near 35 *N. 
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The 10-mb (3l km) temperature (Table ) Increased 10 degrees 
from December 24 to December 25 and 5 degrees from December 25 to December 
26, reaching a maximum of 263 IC on December 26. On December 23, there were 
maxima in x-rays, protons, and AE. All maxima except the x-rays 
evidently originated from the solar flare (importance 2 ) that occurred on 
21 December with a daily flare Index maximum of 142. In this case, the 
heating could have been caused by x-rays, or the electrojet, or a 
combination of these. Table does not indicate an especially strong flux 
of x-rays on December 23. On the other hand, the AE index was above normal 
75 percent of the time on Deeenber 23 (Table 6 ),. about two days before 
the marked 10 -mb temperature change of 10 degrees. 



The Stratospheric Warming of December 197Q-January I 97 I 

^ \ 

This major midwinter warming is a typical example of the so-called 

"European-type" warming. For the first time such an event could be followed 
at regular Intervals at high altitudes up to 42 km with available Nimbus 4 
data. In April I 97 O, Nimbus 4 carried into space, in addition to SIRS B, 
the Selective Chopper Radiometer, which measures the thermal emission in 
selected channels (Houghton and Smith, I 97 O). The radiometer has six 
channels, each accepting infrared radiation emitted by carbon dioxide in 
the' atmosphere below the satellite at wavelengths of approximately I 5 urn. 

The radjAtion measured by channel A originates in a layer approximately 
20 km thick centered at the 2 mb (42 km) level, while that measured by 
channel B originates in a layer of similar thickness centered at 20 mb 
(26 km) . Pigs. iSa and c are northern hemisphere maps of channel A radiance, 
while Figs. I 8 b. and d are northern heinisphere maps of channel B radiance 
^ainett et al , 197^) « Table A gives the black body temperature equi- 
valent to infrared radiance for channels A and B. The interpretation in 
t- 2 nns of atmospheric temperature is explained in Ellis et al.,(l97C); as 
a rough guide it may be considered that these ten^eratures are those occurr- 
ing in the region of the 2 mb and 20 mb levels^respectlvely. 

m ^ ' 

Figs. iBa and b show the slWation on December 30, I 97 O. On, channel 
A (Flg.lSa), a warm region with a maximum radiance of 95mW/m^-sr-cm'^ 
was developing at 45°N 65 °E; 3 days earlier this warm region had been at 

25°N 20°E with a maximum radiance of 85 mW/m^-sr-cm"\ A correspond- 

ing warm region was present on channel B (Fig. 18b) with a maximum at 35°N 
50 °E, implying a northward tilt of the warm anomaly,' 
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Table A 

/ 

Equivalent Temperature as a Function of 
Radiance for a Black Body 


Nimbus ^ Selective 
2 **1 

Radiance mW/m -sr-cm 

30 

40 

50 

6o 

70 

80 

^90 

100 

110 

120 ' 


Chopper Radiometer 


Equivalent temperature 

£ 01.0 

213.7 

224.7 

234.5 

243.6 

251.9 

259.8 

267.2 

274.2 

281. 0 



^ Januaiy 4/ 1971 , a great change had occurred in the temperature 
distributions observed by both channels. The warm area described above 
had. intensified, and at the 2 mb level (Fig. I 8 c) had moved north to 70*^N 
90°E with a radiance at the center of 120 mW/m^-er-cm"^. At this time 

the difference in equivalent temperature near the 2 mb level between the 
hottest and coldest region was more than 60 K. The whole northern hemisphere 


was dominated by a wave number one pattern with little contribution from the 
higher wave numbers apparent on December 30. The pattern was similar at 


the 20 mb level, but the amplitude of wave number one was less and the phase 
is about 25*^ farther eastward. This implies a westward tilt of the warm 
anomaly of 1.5° longitude per km height. Subsequently, the warm region in- 
creased, little in intensity, but moved westward. On January 9, the warm 
region had moved to 75 ‘^W 40°E with a westward slope of 3° longitude/km; 
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thereafter the warm region at the 2 mb level moved farther westward and 
decayed in Intensity. 

On channel A on December 30 (Fig. l 8 a)^ there was a warm region at 65 '^W 
lao^^E and a corresponding warm region on channel B (Fig. l 6 b) at 65°W 17 o“E; 
this was what remained of a warming which bad a maximum intensity on Decem- 
ber 26 at 70 N 125 E at the 2 mb level. On December 30 this warm anomaly had 
a Westward tilt of 3 longitude/km while the cold anomaly over North America 
had a westward tilt of approximately longltude/km. By January 2, 

this warm anomaly had disappeared. 

The increase of amplitude of the warming and the westward phase tilt 
with height are to be expected if temperature variations are caused by up- 
ward propagating planetary waves (Dickinson, I 968 ). 

The Nimbus SIRS instrument was described in the section containing 
the Stratospheric Warming of December 1969 -Januaiy I 970 . Labltske (1972) 
uses both the SCR and SIRS data to describe the warming of December I 97 O- 
January 197^* The set of maps In Pig. 19 shows the intense development 
over Asia on 6 and 7 January I 97 I. In Fig. 19a, the 30-mb temperatures 
are compared with the radiances of channel 8 of SIRS. The northwestward 
displacement of highest radiances from the warm region of the Aleutian high 
is very pronounced Indicating a strong development in the upper stratosphere, 
since conditions in the lower stratosphere (lOO and 50 n4>) are known to be 
similar to those at 30 mb. The lowest radiances are also found somewhat 
west of the cold pole. The pronounced wavenumber 1 should be noted. 



On 6 January (Fig. the highest radiances of channel A have 

reached polar latitudes and are found west of the center of VEUcmest. air 
at the 4-5-kni levels though this, of course, can he stated only tentatively 
because of the few rocket observations. Comparing the position of the lowest 
temperatures with the lowest radiances, we note that the dlsti-lbution is not 
the same as in early winter and that it differs from that of the lower strato- 
sphere in Fig. 19a. The center of lowest radiances in Fig. 19b is found east 
of the lowest temperatures which reflects an eastward tilt with height of the 
axis of the cold pole. This agrees with the temperatures at 6o km (Fig, 19c)^ 
where one part of the cold pole is found over the Barents Sea. 

The axis of warmest air slopes with height across the polar region and 
at the 6o km level (Fig. 19c )„ it is found over Canada, right above the cold 
pole of the upper stratosphere. Cooling took place at 60 km over Asia, 

Europe and Alaska> i.e. , above the intense warming area where the stratopause 
descended to about 40 km; and the northerly winds over Canada and Alaska 
indicate that the Aleutian high moved toward northern Siberia, dominating 
the^ circulation over the whole arctic, while the center of the polar vortex 
was displaced toward Iceland. 

Fig.^0 is a plot (solid line) of the daily maximum temperature at 10 nib 
(31 km) occurring at the higher latitudes in the Northern Hemisphere from 
15 December 1970 to 24 January 1971^ derived from NWP maps . The dashed line 
connects the dally zonal means of channel A of the SCK on Nimbus 4 converted 
into equivalent blackbody temperatures for an average value for and 80 ^ 

N. latitude. (Labitzke and Barnett, 1973) This represents the best estimate 
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available of the high latitude temperature at 2 mb (^3 km). From Fig. 20, 
the 2-mb temperature is greater than the 10-mb temperature from 22 Dec. TO 
to 24 Jan. Tl. ^he 2-mb temperatui'e rises continuously from IT December to 27 
December* Both the 10- and 2 -mb temperatures show maxima near 28 Dec. TO 
and 10 Jan. Tl* The 10-mb curve indicates strong maxima at 4, 8, and 11 
Jan. Tl; the smoother 2-mb curve does not indicate these peaks. However, 
this is not surprising in view of the type of averaging employed (zonal 
averaging). Evidently, the start of the warming is I7 December. 

Fig. 21 shcws.the location.of the warming centers -(maxinrum temperatures) 
at 3-day Intervals on the lO-mb (31 km) WWP Unit charts for 1-13 January 
1971* The crosses indicate the location; the dates are beneath the crosses. 
Two main centers are apparent on January 1 at widely separated' locations; one 
center near 40° K, the other near 60®N. OThese centers merge near 80*N by 13 
January. 

The December 1970 - January 1971 warming represents the first occasion 
in which it is possible to determine in sortie detail the three dimensional 
structure of a stratospheric warming up to 60 km. The period selected 
(Jan. 6-7, 1971) is near the middle of the major warming interval '(Fig. 20). 
Northern hemisphere isotherm patterns were- derived at every 5 km from 25 
to 60 km by using the maps of Fig. 19 (Labitzke,. 1972 ) and the temperature ' 
profiles; in Fig. 22 for, January at 30“N, 50*N, and T0®K (Groves, I9 TO). 

The dashed curve in Fig.22 represents the warming profile for. January 6, 

1971 at 60 derived' from the maps in Fig. I9 and the 10-mb NWP Unit chart. 
Hence, the 25-60 km maps were derived by graphical addition and subtraction. 
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using Fig. 19 and extrapolation, using Fig. 22 , the solid curves of Fig. 

22 for the cold low areas and' the dashed curve for the warming areas. 

Fig. 23; shows the location of the centers of the warm- areas (crosses) 
and cold areas (circles) from 25-6o km on 6 Jan. 71* The altitude in 
kilometers is indicated below the circles and crosses. The closed cold 
isotherm center appears at all altitudes from 25-60 km, while the warm 
center appears only between 30-50 km. The cold center has a decided slope 
with increasing altitude, while the warm center has only a slight slope. 

In addition to these two features, there was another warm center from 
35-55 km located near 25 *N, 240 *E, which was about 5 degrees warmer at 
the same altitude, from 40-50 km, thaathe higher -latitude warm center (Table 

T)*’ The 10-mb NWP Unit charts were not used in this analysis in order 
to maintain internal consistency. The two triangles on Fig. 23 indicate 
the location of two warm centers appearing on the 10 -mb (31 km) chart; the" 
square indicates the 10 -mb low center position. 

Although the information in Table Y is useful, perhaps a more meaning- 
ful way of depicting the December 1970-January 1971 stratospheric warming 
is 'to find those areas of the atmosphere where the temperature differed 
from normal and the magnitude of the difference. In order to accomplish 
this, the' "normal'’ situation must be defined first. To do this, mean 
monthly 10-mb and 5-mb maps for January 1966 (Labitzke, I 968 ) were used 
inasmuch as no major warming occurred during this month. At 60 km, the' 
map for 27 January I 97 I in Fig. 10c of Labitzke (.1972) was used since the 
warming has disappeared by this date and the temperature field has returned 
to normal. The maps at every 5 kilometers from 25“60 km were derived by 
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graphical addition^ subtraction and extrapolation in the same way as 
the "warming" maps above. Finally, temperature change maps were derived 
at every 5 kilometers from 25-60 km and 5 degree isotherm intervals by 
graphically subtracting at each altitude the "normal" map from the 
"warming" map. Fig. 24 shows the temperature change map at 25 km for _ 

6 January 19T1* Table 8 displays the location above K of the tempera- 
ture change centers at every 5 km from 2^-60 km. The maximum increase 
(+45°) occurs from 40-50 km near 80° N between 70-125^ E. There are two 
temperature increase centers at 35 hm and 55 km. There is only one tempera- 
ture decrease center above 40° N-; it exists- between 50-6.0 km*.,. ^.5 

Indicates the location of the centers on a Northern Hemisphere map, where 
a cross denotes a center of temperature increase and a circle is a center 
of temperature decrease. A comparison with Fig. 23 points out that a 
maximum temperature increase or decrease is not located at the same place 
as a warm or cold center. However, evidently there is some correlation;- 
for example, the movement of the temperature decrease center from 50-60 
km (Pig. 25 ) appears to be related to the movement of the cold center 
(Fig. 23 )• Pig* 26 is a plot illustrating how the maximum + AT changes 
with altitude at the same location (78° N 70'^ E; Table 8). The data points, 
selected from the AT charts from 25“60 km, are fitted with two straight 
lines meeting at - 45 km. Fig. 27 shows how the gradient of- + AT varies 
from the center at 78° N 70 ° E along the 70 ° E meridian. The gradient 
can be represented by two straight lines fitted to data points selected , 
from the 45 km AT chart. One line indicates a constant A.Tof 45°K from 
78° N to 66° -N, The other line shows that AT decreases linearly to 28° N. 
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Correlation of Geophysical Indices. December 1970-Janua.ry I 971 

mghllghts for 05-1^ December 1970 : Nineteen events of class-. M or - ■ 

greater were observed during the past seven days - fourteen small class M 
(class M3 or less), one each of class M5, M 7 , and M 8 , and two class XL. 

In addition, a minor proton enhancement and a major geomagnetic storm were 
observed. McMath regions 073 (N15W5^, L :=034) and O 77 (KIOWI 7 , L = 357 ^ 
includes region O 79 which was bombined with 077 on 14 December) accounted 
for 15 of the events. Region 077 produced the following major events - class 
M5 imp in at 09 / 0758 Z with minor radio frequency emission and significant 
ionospheric disturbances; class XI imp IB at II/IOI 6 Z with a 630 flux unit 
burst at 2700 MHa and significant SID; class XI imp SB at I 2 / 0856 Z with minor 
radio emission and major SID; and class M 7 imp Sn at 12/2348Z with minor 
radio emission and moderate SID. Region 073 produced the class M 8 imp 2N 
event near KI 6 WOO at 11/2215Z. 


Hishlights for the period 16-22 December IQ 7 O : Only three events of 


M intensity have 

been reported 

the past seven days: 

I 8 / 0 II 8 Z 

Ml 

No optical counterpart 

22 / 124 5Z 

M 6 

Optically IB, W22W44 

22 / 1516 Z 

M5 

OpticaUy IB, W24W44. Accompanying 

10 cm burst of the order of 100 flux units. 


In addition, a major radio burst with no optical counterpart was observed by 
three different observatories at I 7 / 0549 Z. The peak flux values at 21 and 
11 cm as observed at Carnarvon, Australia were 249 and 5300 flux units, 
respectively. Moreover, Solrad-9 observed no x-ray burst corresponding to 
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this radio burst. Both energetic flares on 22 December arose within McMath 
region 084. The associated spot group made its Initial appearance on_ 18 
December near W24, L = 290 and grew rapidly to a type/area D/18O by I9 
December. Another spot group within McMath plage 084 attained a elass/area 
of D/490 but produced no energetic events. Proton counting rates near earth 
have remained at background. 

It was shown (Fig, 20 ) that the 2-mh zonal temperature rose contin- 
uously from 17-27 December. From this, it was estimated that the warming 
started about I7 December I970. Table 9 presents the indices vs. date 
from 8-24 December 1970. The solar proton maxima are due to electrons 
from the earth^ s radiation belt,. The only solar flare of importance 
2 or'greater during this period was on 11 December, 2215 UT. This flare \m,s 
accompanied by a major radio burst of 1020 flux units at 27OO MHz and 
moderate SID; The ATS-1 satellite observed a proton enhancement with a 
gradual beginning at about I2/050OZ. The 21-70 Mev energy range peaked at 
two particles per sq cm per sec (about four times background) near I2/170OZ. 
Peaks occurred in the 5-21 Mev flux (250). Nighttime absorption of 0.6 db , 
was observed on the Thule 30 MHz riometer. The 14/0200Z peak in the lower 
energy particle flux had a very sharp onset and rapid decay and coincided with 
a geomagnetic sudden commencement at I4/0156Z. 

fl]n Table 9> on 12 December there were maxlJna in the daily flare index, 
solar x-rays, and 280O MHz solar flux. On 14 December, a magnetic storm 
occurred. On this day, there were very pronounced maxima in A , Dst, and AE, 





If the start of the -warming wsase associated with the flare of 11 December, 
then the large AE maximum of 416 on 14 December indicates that the most 
likely avenue was Joule heating associated -with the auroral electrojet near 
100 km. The maximum hour value of the AE index (527) occurred on 14 Decem- 
ber at 2300 UT. The AE index was above the December I97O monthly mean of 134, 
91 percent of the hours on 14 December and 75 percent of the hours on 15 
December. Thus, if electrojet heating were the cause, it took about t-wo days 
to reach the 2-mb surface near 43 km. The start of the warming at 10-mb 
is estimated to be I8 December. The daily 10-mb WWP charts for December I97O 
plainly indicate that a small 243°K isotherm center located at 42N 215E on 
18 December expanded greatly as it moved North to 6o®N, so that it covered 
most of Siberia on 27 December. 

Highlights for the period 23*29 Decen&er ,1970 :' Solar activity was 
at an extremely low level throughout the past seven days. One class M5 west 
limb event was observed from McMath region 084 (Wl8 L = 290) at 26/0840Z. 
Although the x-ray burst was of M5 intensity Its total duration was only about 
10 minutes. The associated radio burst was about 75 flux units at 2700 MHz. 

A minor sudden ionospheric disturbance was reported. A small satellite level 
proton enhancement began at 24/07^0Z, The peak fluxes observed by ATS-1 
occurred at 24/2200Z and were of the order of 40 (5-21 Mev) and eight (2I-7O 
Mev) particles per sq cm per sec. The specific cause of this proton enhance- 
ment is unknown; the nearest energetic events were the class M9 (22/1516Z) 
and M6 (22/1245Z) events near N22"W44 in region 084. The particle enhancement 
may have resulted from coronal storage of particles energized by the above 
events or from an event behind the west limb. The fluxes are still very 
slightly enhanced as of 29/17OOZ, The geomagnetic field was relatively 
quiet throughout the period. 
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Highlights for the period 30 December 197 Q -05 January 19T1 ' No 
energetic flares were observed during the 7“day period. On the average 6 
email regions were visible each day; no one disk feature exceeded 200 
m1,l,Uattthfl In area. A sector-like magnetic disturbance began about 02/0^i-45Z 
followed by a 120 gamma bay roughly 4 hours later. Ho sudden conunencement 
marked the storm’ s onset, however, the co-rotating proton stream was detected 
by ATS-1. 

Highlights for the period 06-12 January 1971 • Solar activity continued 
at a very low level with an average of only ten subflares per day. The most 
energetic event was a class Ml subnormal flare at IO/I 728 Z in McMath region 
111 (S04W32, L=s003). This region dominated the solar disk throughout the week 
with a large (e/ 800 millionths) but magnetically simple sunspot group. It 
was probably the source of a 20-200 MHz radio noise storm which has been in 
progress since O 7 January. The geomagnetic field was also very q.uiet. 

Table 10 lists the indices vs. day of the month from 25 December 197^ 
to 10 January 1971. The E-W and N-S average dally wind components in the 
meteor region above College, Alaska (Hook, 1972) are presented at the bottom 
of the table. A significant change in prevailing winds at the 75-100 km level 
at College was observed, by Hook during this warming. These winds for the 
upper mesosphere, and lower thermosphere were derived from radar measurements 
of the drift of meteor trails. The normal, prevailing wind at College in the 
wind is eastward and usually poleward. During the stratospheric warming, 
both the meridional and zonal prevailing flow reversed from about December 
3l to January 1 , as indicated by the minus values in Table 10. A reversal 
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in the wind was not evident at 130 km, Hencej the upper limit of the circu- 
lation breakdown does not appeeir to have extended above approximately 
100 km at College for this warming event. This wind data by Hook is impor- 
tant because it indicates for the first time the vertical extent of a fully 
developed stratospheric warming. The highest map of Figure 19 goes up to 
60 km; the Hook wind data extends the warming effect upward another 40 km. 

Hook { 1972 ) defines the warming by means of a Temperature Index 
Parameter (TIP)^ which indicates the mean temperature gradient from the pole 
to 60 % and a circulation index^ called the Height Index Parameter (HIP), 
at (Johnson et al. , 1969). Hook concludes that the warming, 

indicated by the minus wind components, started approximately 10 days 
before the wanning at the lO-mb level,, indicated by plus values of TIP 
starting about 10 January. Prom this, he states that the circulation break- 
down during this warming affected the meteor region before it affected the 
mid-stratospheric level. However, the daily 10-mb maximum temperatures In 
Table 10 indicates a minor peak in the warming on 27-28 December. The 
2-mb zonal temperature has maxima on these same days. Based on this, some 
warming apparently occurred before the wind reversal pointed out by Hook. 

Thus, Hook's conclusion that the warming affected the meteor region before 
the mid- stratosphere is questionable. Also, the value of TIP and CIP is 
questionable for indicating the start of a stratospheric warming. As 
pointed out previously, the daily lO-mb NWP charts for December I970 plainly 
indicate that a small S43°K isotherm center located at 42N 225E on I 8 December 
expanded greatly as it moved north to 60 W so that it covered most of Siberia 
on 27 December. It is not easy to pinpoint the start of the warming, but 
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It was indicated previously from 2-mb data that It started near I 7 December 
(Figure 20) . 

Wo flare of importance 2 or greater occurred during the period covered 
in Table 10. The dally flare indices were all below 85 . The low 
level of solar activity is also exhibited in the small solar x-ray flux. The 
2800 MHk solar flux, representing the. EUV, varied in a fairly regular manner. 
The solar proton maxima occur at regular intervals, indicating the influence 
of the trapping region. The maxima in A^, Dst, and AE on December 28 seem 
to have little influence on the 2 -mb and 10 -mb temperature maxima of the same 
date, inasmuch ais these latter maxima did not change from'the previous' day* 
Both temperatures decrease from 28 to 29 December. Table 10 Indicates that 
nothing occurred on 25 or 26 December that could result in the temperature 
maxima on 27 - 28 ' December. Table 9 shows on 22 December an SID with impor- 
tance of 2. and maxima In the dally flare index and solar x-rays. On 24 
December, a ''dlsturbed”day, a small solar proton enhancement was observed. 
Also, maxima occur in A^, Dst, and AE. Since the 2-mb temperature exhibited 
a 3-degree increase from 25 to 26 December, and the. 10-mb temperature showed 
a 5-degree Increase from. 26 to 27 December, electrojet heating is not ruled 
out as the cause for the minor warming peak on 27-28 December. 

Table AO (including Hook* s wind data) show that the stratospheric 
warming Involved the whole atoosphere from about 3o km ( 10-mb) to 100 km 
from about 3l December I 97 O to 9 January I97I. Although it is not possible 
to seek a trigger mechanism during this period, it is possible to look for 
enhancements. For example, the 10-mb charts show maximum temperatures on 4, 
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&nd XX January X9YX. On ths o'th©!' liand^ "the 2-jnb zonaX ■tempera'tures 
Indicate a, fdirXy sbeady rise from 243 K on 3X Deceniber "to a maximum of 
268°K on 10 Januaiy, Therefore, no clues based on enchancement are to 
be sought here. Unfortunately, the AE index was not available for January 
1971, so no check could be made on the possible effect of Joule healiin^;. 



2. CORREILATIONS BETWEEN GEOPHYSICAL INDICES AND STRATOSPHERIC WARMINGS 

A statistical analysis was made to compare average values of the 
indices for warming and nonwarming days with that for total days considered. 
Also., comparisons were made with average values of the indices for 0 to 4 
days preceding the date of start of warmings, and date of major and minor 
warming peaXs. The daily flare index was not Included because it was not 
available in April - August I969. Also, the 2800 MHz solar flux index was 
left out, since the He II index is a more accurate index for EUV. The 
following dates were selected as the start and end of the warmings I969-I97I: 

15 April - 10 May 19^9 
3 July - 21 July I969 

I August - 11 August 1969 

II August - 19 August 1969 

17 December I969 - 4 January 1970 

17 December I97O - 21 January 197I, 

The AE Index 

AE indices were selected during the above warming months. There 
were 98 "warming" days and l34 "nonwarming” days. AE averages for both 
warming and nonwarming days did not differ significantly. Table 11 dis- 
plays the AE index versus the elapsed days from start of the warmings above, 
for up to 4 days preceding the warming. The April - May I969 warming was 
not Included Inasmuch as the start time could not be determined accurately. 
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In all five cases, AE had a large maximum one to three days preceding 
the start of the warming. Table 12 shows the relative variation of AE 
from the average value of each 5-day period vs. elapsed days from start of 
warming. There is marked variation in the day -by-day values of AE in the 
days immediately preceding the start of the warmings. The average value of 
AE for the two days preceding the start of all warmings was 46 percent above 
the averages for the five-day Intervals considered; it was 37 percent below 
average on the day the warmings started. The relative AE variations are 
even greater, for all individual warmings, except for that starting 11 August 
1969 . The others show variations ranging from 59 percent to 183 percent 
a'jove average at 1-3 days preceding the start of the warmings. It is not 
unusual to find variations of the type indicated in Table 12 during any 
five-day Interval. However, it may be significant that the daily variations 
line up in the fashion described, where on 0 day, all AE Indices are below 
normal, the -.1 day, one is below normal, the -2 day, none, -3 day, three are 
below normal, and the -4 day, four are below normal. 

The relationship of the AE index to warming was explored further with 
regard to the major and minor warming peaks. Table 13 exhibits the AE 
index vs. elapsed days from date of ten major warming peaks. The AE rela- 
tive maxima are marked with an asterisk. It is noted that there are twice 
as many relative maxima on either -land -2 as on 0 and -3 days. The 
average AE value Is a maximum at -2 days (222). The Grand AE average (I 98 ) 
is 14 percent above normal. Table 14 shows the AE index vs. elapsed days 
from date of 12 minor warming peaks. Again, as in Table l3 , the maximum 
AE average occurs at -2 days. Over twice as many relative maxima in AE occur 
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on -2 day ae on either 0 or -4 days. The average AS (149) for 0-4 days 
preceding minor warming peaks Is below normal. In Table l3, the average 
AE for 0-4 days (total 50 days) preceding major peaks (I 98 ) is one-third 
higher than the same average for the minor peaks (total 60 days). This 
may be significant since it was found previously that there was no signifi- 
cant different between the AE average for warming days versus nonwaimiing 
days. 

Solar X-Rays 

The average msocimum hourly value (24 hour interval) of the 1-8A 

_3 - 2-1 • 

solar x-ray flux is 7»9 ^ 10 ergs cm s . The average value is less 
during the II 9 warming days than during the 144 nonwarming days- Table 15 
shows the x-ray flux values versus elapsed days from start of warming- 
The average flux occurring during days 0 to -4 is less than average. The 
average flux is greater than average on day 0; also, there are three 
relative maxima. Tables ^ 3.6 and' I 7 exhibit the solar x-ray flux vs. elapsed 
days from date of major and minor warming peaks, respectively. In both 
tables, the average x-ray flux occxirring 0 to 4 days preceding major and 
minor peaks is below average. Nothing unusual is noted in the distribution 
of the relative maxima. 

Solar Flares 

Forty percent of the flares of importance 2 or greater out of a total 
of 69 occurred during warming days. Since 45 percent of the total days 
considered (total 264) were arming days, apparently, there is no significant 
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difference in the occurrence of flares (importance s 2) hetween warming 
and nonwarraing days. Only 2 out of 10 importance 3 flares occurred during 
warming days. 

Table 18 displays the solar flare importance (2 or 3) versus elapsed 
days from start of wanning. The number of flares occurring 0 to 5 
d^s preceding the start of the warming Is not significantly different 
from the average. However, in 2 out of 5 cases, flares of importance 3 
occurred at -3 days. Tables 19 and 20 show the solar flare importance 
versus elapsed days from dates of ma^or and minor warming peaks, respec- 
tively. In both tables, the number of flares occurring 0 to 5 days preceding 
major and minor peaks are not significantly different from the average. 

Nothing unusual is noted in the flare distribution, 

Tbe Index . 

The average value of for the interval considered (264 days) was 10.0. 

The average (9*0) for II 9 warming days is less than the average (10.8) for 

the 145 nonwarming days. Table £1 shows the A^. indices vs. elapsed days 

from start of warming. Above average values of occur at -2 and -3 days. 

All the relative maxima of A occur at -1 to -3 days. The average A for 0 

P 

to 4 days preceding the start of the warmings is below normal.. Tables 22 
and 23 exhibit A^ vs, elapsed days from date of major and minor warming 
peaks, respectively. The average for 0 to 4 days preceding the peaks is 
below normal in both tables. Both tables show a relative mximum at -2 days. 
This maximum is above normal (28 percent) for major peaks and below normal 
for minor peaks. 
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The Index 

Figure 3 (Part I) shows the relationship between and AE. The 
main feature to be noted is that the main phase of a magnetic storm has 
high positive AE values but large negative values. Hence, in the 
following analysis, " m a x i m a” in are taken to be the largest negative 
values. The average index (26l days) is -0.6. The average for II9 
warming days (+O.7) is below normal; the average for 142 nonwarming days 
(-2.2) is above normal. Table 24 displays the D . index vs. elapsed days 
from start of warming. is above normal at -1 to -3 days. The relative 

V '-maximum" is at -I days. The, average from 0 to 4 days preceding the start 
of warming is below normal in 4 out 5 cases, tables 25 and 26 show the elapsed' 
days from date of major and minor peaks, respectively. ' Both tables show 
that the average 0 to 4, days preceding the peaks is below normal. The 
relative maximum in Table 25, major peaks, occurs at -2 days, while It occurs 
at -3 days, Table 26, minor peaks. 

The He II Index 

This index, representing the 304A solar flux, is a more accurate 

indication of the EUV upper-atmosphere heating than the 2800 MHz solar 

flux (Anderson, 1965)* However, it is available only for the months of 

April - August 1969. During this period, the He II index ranged from 68 to 
-2 -1 8 

96 photons cm s x 10 . There was no significant difference between the 
average for the 64 warm days (85»8) and the average for the 75 nonwarra 
days (84.3). The average for all days considered was 64.9-. In -Table 27, . 

He II index vs. elapsed days from start of warming, the relative inaximtun 
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occurs at 0 days. The average He II index for 0 to 4 days preceding the 
start of the warming is 8o, a helow-normal value. Tables 28 arid 29 'exhibit 
the He H index vs. elapsed days from date of major and minor warming peaks, 
respectively^ In both tables, the average of the He II index for 0 to 4 
days preceding the peaks is close to normal. 

Solar Protons (greater thaji 10 Mev) - ^ 

Every three to four days a large maxinrum appears in the solar proton 
flux, evidently due to passage of the measuring spacecraft through the radia 
tion belt (see, for example, Table ^lo). Therefore, it is meaningless to 
attempt to analyze the proton flux like, say, the x-ray flux (Table 15 ). 
However, a check was made to try to find afly influence, if any, of the 
radiation belt particles on stratospheric warmings (start, n»Jor' and minor 
peaks). Hone was found. 



The SID Index 


An equal number (24) of vanning and nonwarrolng days were associated 
with an SID index of 2- or greater. Nineteen percent of the days had an SID 
2- or greater. Table 3o shows the SID Index vs. elapsed days from start of 
warming. Twenty percent of the days 0 to 4 days preceding the start of 
warming had an SID 2- or greater; hence j there is nothing' unusual about 
this distribution. 

Tables 3l and 32 display the SID index vs, elapsed days from date 
of major and minor warming peaks, respectively. In Table 31, only 9 per- 
cent of the days 0 to 4 days preceding the date of major warming peaks have 
an SID index of 2- or greater. In Table 32, 28 percent of the days have 
such an index. Consequently, about three times as many SID indices 2- or 
greater are found in the 0 to 4 days preceding minor warming peaks than in 
the same period for major warming peaks. 

■Magnetic Condition 

A magnetically-disturbed day is denoted by the indices D (disturbed) 
or S, (principal magnetic storm). Nineteen percent of the warming days had 
a D or S index, while 19 percent of the nonwarming days had such an index. 
Nineteen percent of the days considered (253) had either a D or S index. 

Table 33 shows the magnetic condition index vs, elapsed days from start of 
warming. The number of disturbed days is about average. Tables 34 and 35 
exhibit the magnetic index vs. elapsed days from dates of major and minor 
warming peaks, respectively. The number of days in both tables having a 
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magnetic index 0 to 4 days preceding a peak is I 7 percent, or almost 
average. In Table 35? 91 percent of the days having an index occur within 
2 days preceding the minor warming peak. 



The preceding analysis indicated no strong correlation of solar flares, 
SID, x-rays, He II and magnetic condition with stratospheric warmings. 
However, some definite correlations with stratospheric warmings were in- 
dicated for the magnetic indices AE, A and D . , These magnetic indices 

p st 

are correlated somewhat with each other; for example, the index 
contains contributions from at least two major sources, the auroral 
electrbjet (AB) and the ring current (D . ), (See Figure 3, Bart I and 
Figure '2a, b, Coffey and Lincoln, 1973.) (These indices are described 
in greater detail In Baxt I. ) 


^ Table 36 summarizes the statistical information concerning AE, A^ 
and It shows the average value of the indices vs. warming, non-warming 

days, total days, and for 0-4 days preceding the start of warming, major 
wanning peaks and minor wanning peaks. The symbols in parentheses after 
the averages refer to whether or not the index average is above-normal 
(A), below-normal (B), or normal (N). 
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In Table 36, there is no significant difference between the AE index 
averse for warming days and nonwarming days. Althovjgh the average AE for 
0 to days preceding the start of warming is 2k percent below normal in all 
5 cases considered, AE has a large maxlinum 1-3 days preceding the start of warming 
ra n g ing from 20 percent to 100 percent above average. The average AE for ail 
cases at 2 days preceding the start of warming la k6 percent above normal, 

^dille at 0 and -U days it is 37 percent below normal and 2k percent below 
norm^, respectively. The average AE for 0 to days preceding 10 major warming 
peaks is 1^ percent above normal. The maximum average AE is at 2 days preceding 
the major peaks, where it is 2? percent above normal. The average AE for 0 to 
k days pjrecedlng 12 minor warming peaks is 15 percent below normal. The maximum 
average AE Is at 2 days preceding the minor peaks. 


Altto\agh the average A^ for 0 to k days preceding the start of warming Is 
13 percent below normal (Table 36), the average values at 2 and 3 days preceding 
the start of warming are 20 percent and TO percent above normal, respectively. 
The average A for warming days is lower than for nonwarming days. The average 

ir 

Ap at 0 and k days are 66 percent below normal and 3k percent below normal, 
respectively. The average for 0 to *1- days preceding 13 major peaks is 9 
percent below liormal. Ttie m a x imum average is at 2 days preceding the major 
peaks, where It Is 28 percent above normal. The average preceding 13 minor 
warming peaks is 27 percent below normal. The maximum average is at 2 days 
preceding the minor peaks. 


11-38 



In Table 36, the average for vajrmliig days is lower than for nonwarming 
days. The average for 0 to 4 d£^ preceding the start of the warming is 
below" noinal for 4 out of the 5 cases considered. The average is "above" 
normal at 1 to 3 days preceding the start of warming, with the maScimum average 
at -1 day. The average for 0 to 4 days preceding both the maior and minor 
warming peaks is below normal. The wia.YtnniTw average value occtirs at -2 days for 
major peaks and -3 days for minor peaks. 

Thxis, there appears to be a significant relationship, between the AE index 

1-3 days preceding the data of the start of the warmings and dates of the major 

and minor warming peaks. The relationship Is strongest for stai^ of the warmings 

and weakest for minor warming peaks. ’She analysis of the A index results in 

P 

the same concltision. The analysis of the index does not give as strong 
results. It is of interest to note that none of the indices in Table 36 had 
above-normal averages for warming days. 



3. DISCUSSION AMD CONCLUSiCWS 

Five different stratospheric warmings, three different associated 
pbc:nomena (start, major and minor peaks) and three different magnetic 
indices all indicate the same thing; namely, the average value of the 
magnetic indices {AE, Ap, display a marked increase 1 to 5 days before 

the warming phenomena and a marked decrease at 0 and -ii days. Although the 
temperature of the thermosphere above 100 km is definitely enhanced during 
magnetic storms (Anderson, 1973)^ do direct evidence has been presented 
that the mesosphere and lower atmosphere is affected. Statistical studies 
have established that correlations exist between lower atmospheric behavior 
and variations in the Intensity of , georaegnetic fluctuations (Shapiro, 1972; 
Roberts and Olson, 1975; Wilcox et al., 1973). King (19?4), based on a 
comparison between meteorological pressures and the strength of the geo- 
magnetic field, suggests a possible influence of the Earth's magnetic' 
field on the height of the 500 mb level in regions surrounding the magnetic 
poles. These statistical studies have indicated that indirect correlations 
exist between magnetic variations and lower atmosphere circulation. 

However no direct (cause and effect) correlations can be established 
unless a physical theory is advanced to explain how the connections work. 
Hines (l973) points out that the observed magnetic fluctuations are un- 
doubtedly to be associated with solar plasm mechanisms. The existence 
of fluctuating winds in the ionospheric dynamo region and of associated 
magnetic variations at ground level is equally beyond question. The main 
question is: are these dynamo winds influenced by the lower atmosphere 

circulation? In other words, does the lower atmo^here circulation cause 


changes in the ionospheric winds so that magnetic variations ensue or 
vice versa? 

The major stratospheric warmings from 1951 to 1966 all occurred 
dxiring years of strong solar and geonmgnetic disturbance, or following 
active bursts of disturbance (Willett, I 968 ), The only six years lacking 
stratospheric warmings were also the six geomagnetically quiet years. 

The definite correlation of the warming phenomena with magnetic activity, 
especially the AE index indicates that the auroral electrojets near 100 km 
may be the primary source of the auroral zone heating that influences the 
circulation at , 50 . km and below. 

The auroral electrojets are intense electric currents, associated 

with polar magnetic disturbances, that flow in the auroral oval (Akasofu 

et al., 1965). The westward auroral eledtrojet grows to approximately 

6 

IQ an^eres during a magnetic substorm due to an increase in the conduc- 
tivity of the E-layer of the ionosphere in the midnight to dawn sector of 
the auroral oval. The eastward flowing electrojet has a strength of about 

10^ anaperes in the late evening sector. Cole ( 1962 ) has estimated that 

- 2-1 

5 ergs cm s could be deposited in the auroral zone by means of joule 
heating associated with the electrojet at moderately disturbed times and 
an order of magnitude more during large geomagnetic storms. Whether or 
not the electrojets are the primary cause or trigger the heating from 
below (Hxrota, I 967 ; Matsuno, 1971; Trenberth, 1972) cannot be determined 
without a detailed energy transfer analysis, which cannot be performed at 
present because of lack of data in the mesosphere. 


ii-ia 



Concluaions 


This detailed study of five different stratospheric warmings 
(1969-1971) has established that the magnetic indices AE, and Dst 
display a marked increase 1 to 5 days before warming phenomena and a 
marked decrease at 0 and -4 days. The definite correlation of the warming 
phenomena with the AE index indicates that the auroral electrojet Joule 
heating near 100 km may be the original source of the heating that in- 
fluences the atmosphere at 50 km and below. If this is true, then whether 
the electrojet heating is, the primary cause or triggers the heating by 
releasing energy in the troposphere cannot be determined without meso- 
spheric data above 50 km. 
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n-it-3 



(1-17 July 1969) 


9 

10 

11 

12 

13 

14 

15 

16 

17 

Units 










Importance 

Importance 

1 

3 * 

1 

7 ^ 

.3 

3 

3 

2 

1 

Max. Hr. Value 
lO'^ergs 










cm s"-*- 

91* 

89 

87 

89* 

89* 

85 

85 

85 

84 

10' photons 
cm"^s“^ 

163 *. 

159 

153. 

147 

150* 

139 

122 

122 

120 

-22 -2 
10 yvjm 

1 

1 12812* 

1 

1 

1 

437* 

1 

1 

Max. Hr, Value 
Protons cm "2 
s'^(>lCMev) 




s 








7 * 

6 

12 

14 * 

13 

5 

10* 

4 

Daily Avg, 

8 

5 * 

14 

7 

2 

0* 

6 

6 

6 

Max. Hr, Value 

299* 

235 

208 

291* 278 

246 

129 

229* 120 

DailyAvg. 

rrnmmnfi 

60.5 

60.5 61.0 

60.7, 

61.4* 

60.3 

60.1 

58.8 58.2 

-2 -1 

Ergs cm s 
ster'^cm"^ 


TABLE 


2 


•1-8A Solar X-Rays (Explorer 37) 



, _2 -1 
(10 ergs cm s ) 

Start 

(UT) 

Peak 

(UT) 

End 

(OT) 

Duration 

(Min) 

1 

96 

1125 

1128 

1135 

10 

1 

130 

1627 

1629 

1636 

9 

2 

89 

0349 

0352 

0400 

11 

2 

44 

0408 

0412 

0416 

8 

3 

1900 

1523 

1524 

1603 

40 

4 

61 

0006 

0014 

0032 

26 

4 

38 

1200 

1206 

1206 

6 

5 

,210 

1229 

1237 

1244 

15 

7 

.53 

1636 

1648 

1648 

12 

8 

.46 . 

2218 

2227 

2234 

16 

10 

38 

0649 

0653 

0658 

9 

12 

86 

1430 

1435 

1517 

47 

13 

38 , 

0341 

0348 

0354 

13 

13 

71 

0400 

0410 

0416 

16 

13 

46 

1513 

1515 

1517 

4 

13 

39 

2229 

2245 

2245 

16 

13 

58 

2320 

2324 

2328 

8 

14 

43 

0244 . 

0246 

0249 

5 

14 

61 

1437 

1437 

1438 

1 

14 

40 

1730 

1736 

1741 

11 

15 

57 

1717 

1720 

1728 

U 

16 

37 

1903 

1910 

1914 

11 


IX^kh 


TABLE 3 

Indices vs. Date (29 July - 14 August 19^) 


Date 


Index. 


29 


30 31 


10 


11 12 


13 14 


Units 


Solar Flare 
SID 

Solar X-rays 
(1-8A) 


2 

2 2 


(9 


m 6 


3* 


I 

6*1 


He II Flux 


71 

70 

8Q 

1 

83 

86 

93* 

88 

88 

V 

88 <0 

2800 MHz 
Solar Flux 

• \ 

133 

'l43 

\ 

162 

I 

171 

187* 

183 

.188* 

182 

167 

159 146^ 

■ \ 

Solar 

Protons 

1 

\ 1 8o35» 

1 

1 

1 , 

1 

6* 

1 

1 

45470* 1 ' 

Magnetic 

Condition 


J 


1 

! 


D 




S 

A 

2 


. 6 

■(3 

5 

14* 

13 

7 

5 

8; 

P 

^st 

AE 

-1 

1^ 

0 

16 

21 

9 

-13* 

. -4 

4 

16 4 

67 

©, 

\ 

102 

93 

294* 

291 

178 : 

.116 

215* 210 

Max. Radiance 
SIRS-A, 

57.5* 57.^ 

s. 

57.0' 

56.9 

58.0 

; 

59.3 

60.5 

61,4 

64.2* 63.4 - 63,_6* 


90* 9C^ 


34* 


8T 




(^x 

6 

© 

' 4 

S 

21* 

Nw/ 

-1 


“15 

-3 

203 


'65 

315* 

61.9 

s 

59.6 

I 

i 

"57.2 

58.0 


Importance 

Importance 

2 2 Max.Hr. Value 

10“-^ergs cm“2 

o 

85 82 lO^photons 

cm"2 s“^ 

121 m 10 ‘^Wn'^MHz*^ 

Max, Hr, Value 
1 2221* Protons cm" ■ 
s“^(>10Mev) 


7 6 Daily Avg, 

-4* . 0 Max. Hr. Value 

174 192* Daily Av. 
Gammas 


Channel 8 


59.1 59-4 Ergs cm s 
•— l/vn"i 


-2 -1 

ster"*^Gm* 


* Relative Maximum 


H 

H 

I 

*r 

VJ 1 



.Date 


TABIiE 4 

Indices vs. Date (11-27 December 19^9) 


11 12 13 14 15 i6 IT 18 19 .20 21 


\ 

Index''. 

V 

11 

12 

13 

Solar Flare 




Daily Flare 

3 

11 

. 70 

Index 



SID 

3' 

2 


Solar X-rays 

' 4 

9* 

2 

(i- 8 a) 



2800 MHz 
Solar Flux 

118 

122 

130 

Solar 

Protons 

1 

1 

4262* 

Magnetic 

Condition 




A 

11* 

5 

3 

P 



8 


-5-^ 

-1 

AE 

140* 

88 

46 

10 -mb Max. 

m 

234 

238 

240* 


22 23 24 25 26 2? 


Units 


© 

% 


0 





2 

.26^ 

I 


\ 



16 36 45 22 49 54 

.+ 


10 * 


2 \ 2 
■ \ 


2 

6 * 


2 

7* 


I 


118 122 130 133 ':^36 137 l'^7 15^ 160 163* 156 153 15^* 151 1‘^J 


\ 


1 1596* 

\ 


I 

1 1705* 1 

i 


5 W 

10 


1 '2631* 1 

\ 

'd 

7 ( n -» 


1 2101* 


Importance 


Importance 


Max.hrly. 

value 10“^ , 

-2 -1 

ergs cm 

^^-22 -2 
10 w m 

MHz-1 


Max.hrly. 


value Protons 
cm'^s" 

(>10 Mev) 


10 * 


Daily Avg. 




\ 


Temp. 


10. 

14 

20 

17 

.19 

5 

• 1 

-1 

-2* 

4 

-5* 

Max.hrly. 

value 

47 

39 

32 

35 

! 

107 

(2^^ 155 

189* 

147 

145 

Daily Avg. , 
Gamnas 

248* 

241 

243 

238' 

243 

243 

248^ 

'258 

263* 

258 " 

»K. 


* Relative maximum 



TABLE 5 

1-8A Solar X-Rays (Explorer 37) 



Flux 

Day(Dec. I969) (lO ^ergs cm 


14 

87 

14 

260 

17 

260 

17 

770 

17 

470 

23 

91 

23 

62 


Start 

(UT) 

Peak 

(ra) ' 

End 

(OT) 

Duration 

(Min) 

0257 . 

0305 

0321 

24 

033,8 

. 0351 

0444 

6 

0038 

0052 

0230 

112 

0658 

) 

0700 

',0726 

28 

! 

0950 

. i 

0956 

1017 

27 

p 64'9 

/"■0654 

] 

. 0750 

61 ' 

0850 

\ 

0900 

r 

0922 

32 





TABLE 6 

AE Index (Gammas) (11-27 December I969) 


Date 

Mean Value 

Maximum Hourly Value 

Percent Hours/7110 

11 Dec. 1969 

140* 

284 

58 

12 

88 

221 

25 

13 

46 

109 

0 

14 

62 

260 

9 

15 

99 

278 

42 

16 

190* 

521 

.. 71. 

17 

, 78 

176 

25 

18 . 

47 

113 

9 

19 

39 

86 

0 

20 

32 

■ 55 

0 

21 

35 

83 

0 

22 . 

107 

364 

42 

23 

236* 

527 

75 

24 

155 

372 

42 . 

25 

189* 

324 

67 

.26 

147 

, 413 ' 

54 

27 

145 

429. 

.58 


* Relative maximum 
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Table 7 


LOCATION AND TEMPBR/OTRE OF NOBTHEEN HEMISPEERE 
COLD AND WARM CENTERS, 25“60 KM, JAN. 6, 1971 



Cold Center 

Warm Centers 



Altiti^de 

(km) 

Lat . 
(®W) 

Long. 

(°E) 

Temp. 

(°K) 

Lat. 

(°N) 

Long. 

(°E) 


Lat', 

Long. 

Temp. 

(°K) 

Lat . 
(°N.) 

Long,' 
(‘^E) ; 

Temp, 

(°K) 

25 

TO 

10 

198. 

60 

150 

233 

30 

15 

218 




30 

TO 

10 

203 

62. 

.155 

2-43 

30 

15 , 

228 




35 

T 8 

320 

228 

58 

145 

253 

30 

15 

243 ' 

25 

240 , 

243 

Ao 

6 o 

240 

233 

65 

l 4 o 

268 




25 

. 240 ' 

273 . 

45 

57 

245 

243 

65 

140 

278 




25 

245 

283 

50 

60 

250 

248 

65 

,l 4 o 

- 288 




25 

240 

293 ^ 

55 

TO 

220 

248 




, 



30 

250 

278 

60 

TQ 

65. ' 

10 

200 

238 , 
233 



' 

55 

250 

278 


■ 
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Table S 


LOCATION AND MAGNITUDE OF NORTHERN HEMISPHERE 
TEMPERATURE CHANGE CENTERS, 25-60 KM, JAN. 6, 19T1 












TABLE 9 

Date (8 - 24 December l970) 


15 

16 

17 

18 ■ 

19 

20 

16 

10 

6 

61* 

52 

16 

2 

2 

1 

4 

2 

2 

■54 

147 

147 

152 

152 

156^< 

4 

1 

.1 

13522* 

1 

1 




S ' ■ 



14 

4 

3 

3 

10* 

8 


50 2 


Importance 


Importance 


2 2 1 Max. Hr. Value 


4 1 .1-13522^^ 1 1 9957* i i Max. Hr. Value 

'V' Protons cm*^ 
s"^C> 10 Mev) 

S' ■ ‘i)' 

14 4- ■ 3 3 10* 8 4 6 7 Daily, avg. 

-56 -35 -9. 3 -17* -14 -9 -2 Max. Hr. Value 


255 -71 ^5^ 47 - 116 186* 104 130* 98 l^; Da 

^ I 

236 236 236 238 , 240 241 243 243 244 245 ^ 

238 243 . 243 ^3. 243 243 243 243 243 243 °K 


^ Daily Avg. 
Gammas 











TABLE 

10 














Indices vs 

, Date (25 December 19?0 

- 10 January 1971y 






^'N^ate 



















Index 

25 

.. 26 

27 

28 

29 

30 

31 

1. 

2 

.3 

■ 4 

■ 5- 

6 

7 

8 

y 

10 

Units 

Solar Flare 


















Importance 

Daily Flare 

85 * 

20* 

1 

0 

4 

12* 

0 

10 

6 

0 

5 

18 * 

10 

10 

1 

65 * 

12 


Index 


















SID 

2 
















2 

Importance 

Solar X-rays 
(1-8A) 

1 , 

3* 

1 

0 

1 

1 

1 

1 

1 

1 

1 ‘ 

1 

2 * 

1 

1 

4 

5^ 

Max. Hr, Value 
10 "^ergs cm'^ 
s“^ 

2800 MHz 
Solar Flux 

124 

122 

120 

117 

123 

129 

133* 

130 

135 

135 

140 

146 

147 

i5C^ 

149 

153 * 152 

10"^^w m“^z'^ 

Solar 

Protons 

8946 * 

' 3 

2 

4242* 

1 

1 

6429 * 

1 

1 

1 

9595* 

1 

1 

13009 * 

1 

1 15600 * 

Max.Hr.Value 


















Protons cm“2 
b“^(> 10 Mev) 


Magnetic 

Condition 



S 

S 

S 



S 

S 

T> 









A 

P 

4 

4 

8 

18 * 

12 

13* 

2 

6 

19 

33* 

16 

9 

4 

1 

0 

1 

6 

Daily avg. 

^ 4- 

st 

-8 

-5 

-2 

-20* 

-17 

-16 

-1 

5 

7 

-8* 

-8*; 

-3 

-2 

3 

9 

17 

19 

Max. Hr. Value 

AE 

71 

74 

179 

324* 

159 

143 

25 











Daily Avg. , 

N-SWlnd, 

College 

+8 

+8 

+8 

+5 

>5 

-20 

. -25 

-40 

-40 

-38 

-38 ; 

-35 

-35 

-35 

+20 

-20 

+10 

Gammas 
Avg. Daily 
C<^ponent , 

(75-105 3on) 

















ms“^ 

2 -mb Zonal 
Temp( 60 - 80 °) 

246 

249 

252 * 

252 * 

250 . 

249 

249 

250 

251 

252 

254 ■■ 

255 

257 

260 

263 

266 

268 * 


10-mb Max. 
Temp. 

243 

243 

248* 

248* 

243 

243 

243 

238 

238 

243 

253 * 

243 

243 

253 

258 * 

248 

248 


E-W Wind, 
College 

+15 

+10 

0 

+25 

+35 ■ 

+5 

-15 

-35 

-35 

-35 

-35' 

-35 

-35 

-38 

-40 

+2 

+15 

Avg. Daily 
Qomponent , 
ms“^ 

( 75-105 Hm) 



















4 * 


Relative maximum 


TABLE 11 


<’ 


AE ; 

INDEX VS. 

ELAPSED DAYS EROM START OF 

WARMING 


Elapsed 

Days 

Start Date 

0 

-1 

-2 

-3 

-4 

Average 

3 Jill 1969 

97 

169 

294* 

125 

» 

126 

162 

1 1969 

102 

131 

191* 

67 

108 

120 

11 Aug 1969 

65 

219 ** 

203 

210 

# 

215 

182 

17 Dec 1969 

78 

190» 

99 

62 

46 

95 

IT Dec 1970 - 

.75: 

71 

255 

. 416 *. 

69 

147 

Average 

83 

156 

208* 

176 

113 

141 


* Relative maximum 
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TABLE 12 


0 

RELATIVE VAIilATION OF AE FROM AVERAGE VALUE OF EACH 5 -DAY 
PERIOD VS. ELAPSED DAYS FRCM START OF WARMING 


Elapsed 

Days 


Start Date 

0 

-1 

-2 

-3 

-4 

3 Jul 1969 

.60 

[ 

1.04 

1.81* 

.77 

.78 

1 Aug 1969 

.85 

1.09 

1.59* 

.56 

..90 

11 Aug 1969 

.36 

1.20* 

. 1.12 

1.15 

1.18 

17 Dec 1969 

.82 

2.00* 

1.04 

.65 

.48 

17 Dec 1970 

.51 

. 48 . 

1.73 

2 * 83 * 

.47 

Average 

.63 

1,16 

1 . 46 * 

1.19 

.76 


* Relative rriaxlinum 


r 







TABLE 13 

A 1 -: VS. ELAPSED DAYS FJiCM DATE OF MAJOR WARMING PEAKS 


Elapsed 


Peak Date 

0 

1 1 < I 

-1 

-2 

-3 

-4 

Average 

30 Apr 1969 

A21* 

207 

333 * 

216 

158 

267 

5 May I969 

235* 

191 

380 * 

342 

222 

274 

8 Jul 1969 

139 

146 * 

98 

110 

90 

116 

11 Jul 1969 

208 

235 

299^ 

139 

146 

103 

13 -Jul 1969 

278 . 

291* 

^ 208 

-235 . 

299 *- 

262 

6 Aug 1969 

116 

178 

291 

294* 

93 

194 

8 Aug 1969 

210 

215* 

116 

178 

291* 

202 

16 ,Aug 1969 

121 

50 

192* 

174 , 

315* 

170 

26 Aug 1969 

147 

189* 

133 

236 * 

107 

167 

29 Aug 1969 

65 

66 

145 

147 

189* 

122 

Average 

19 ^*. 

177 

222* 

207 

191 

198 

* Relative maximum 
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^AELE 14 

A]'; VS. E-LAPSED DAYS FROM DATE OF MINOR WARMING PEAKS 


Elapsed 


Peak Date 

0 

-1 

*2 

-3 

-4 

Average 

19 Apr 1969 

124 

271 

464 * 

363 

319 

308 

23 Apr 1969 

100 

115 

88 

174* 

124 

120 

12 May I969 

157 

157 

211 

218* 

134 

175 

3 Jul 1969 

97 

169 

294* 

125 

126 

162 

18 Jul 1969 

132 

■ 120 

■ 229 ■ 

■■ 129 ' ' 

246 * 

171 

22 Jul 1969 

149* 

138 

121 

78 ■ 

132 

12K. 

20 Aug 1969 

205 

346 * 

115 

230 

121 

203 

13 Dec 1969 

46 

. 88 

140 * 

114 

99 

97 

18 Dec 1969 

47 

78 

190* 

99 

62 

95 ' 

20 Dec 1969 

32 

39 

47 

78 

190* 

7t 

2 Jan 1970 

403 * 

120 

65 

44 

65 

139 

27 Dec 1970 

179 

74 

71 

156* 

98 

116 

Average 

139 

143 

170* 

151 

143 

' 149' 


* Relative maximum 
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TABLE 15 

Solar X-Rays. ( 1 - 8 a) (Max. Hr. Value, 10 *^ ergs cm“^ s"^) 
Vs. Elapsed Days From Start of Warming 


Elapsed 


Days 


Start Date 

0 

-1 

-2 

-3 

-4 

Averag 

3 Jul 1969 

3 

3 

2 

4* 

1 

2.6 

1 Aug 1969 

7 * 

2 

2 

4* 

2 

3.4 

11 Aug. 1969 

10* 

2 

1 

1 

6* 

4.0 

IT Dec 1969 

. 26* 

2 

3 

9 ^ 

2 

8 . 4 * 

17 Dec 1970 

1 

2 

2 

5 

5 

3.0 

Average 

9.4* 

2.2 

2.0 

4,6 

3.2 

4.3 


* Relative maximum 


TABLE 


16 


-Solar X-Bajrs (I^A) (Max. Hr. Value, 10 “\rgs cm'^ s“^) 
Vs. Elapsed Days From Date of Major Warming Peaks 

Elapsed 


Days 


Peak Date 

0 




-4 



-2 

-3 

Average 

30 Apr 1969 

2 

1 

1 

4 * 

2 

2.0 

5 May I969 

1 * 

2 

2 

7 * 

3 

4 . 2 . 

8 Jul 1969 

2 

1 

1 

6* 

4 

2.8 

11 Jul 1969 

1 

3 * 

1 

2 

1 

1.6 

13 Jul 1969 

3 

7* 

1 

3 

1 

3.0 

6 Aug 1969 

1 

3 * 

2 

3 * 

6 

3.0 

8 Aug 1969 

1 

6* 

1 

3 

2 

2.6 

16 Aug 1969 

0 

1 

2 

2 

2 

1.4 

26 Aug 1969 

7 

5 

5 

8* 

6 

6.0* 

29 Aug 1969 

2 

3 

2 

7 * 

5 

3.8 

4 Jan 1971 

1 

1 

1 

1 

1 

1.0 

8 Jan 1971 

1 

1 

2* 

1 

1 

1.2 

11 Jan 1971 

1 

5 * 

4 

1 

1 

2.4 

Average 

' i 

2.2 

2.9 

1.9 

3,7* 

2.7 

2.7 


*. Relative may itniipi 


n -$8 



TABLE 17 

Solar X-Rays (1-8a) (Max. Hr, Value, 10“^ ergs cm"^ s"^) 
Vs, Elapsed Days From Date of Minor Wafmlng Peaks 


Elapsed 

Days 


Peak Date 

0 

-1 

-2 

-3 

-4 

Average 

19 Apr 1969 

5 

4 

1 

9 

10* 

5.8 

23 Apr 1969 

9 

4 

20»^ 

5 

5 

8.8 

12 May I969 

8* 

1 

i 

1 

1 

2.4 

.3 Jul 1969, , 

3 

3 ; 

2 

4* 

1 

2.6 

18 Jul 1969 , 

0 

1 

2 

3. 

3 . 

1.8 

22 Jul 1969 

1 

1 

' 1 ■ ■ 

1 

0 

0.8 

20 Aug 1969 . 

6* 

0 

2 

0 

0 

1.6 

13. Dec. 1969 

'2 

9# 

4 

1 

1- 

3.4 

18 Dec 1969 

2 

26 ^ 

2 

3 

9 

8.4, 

20 Dec 1969 

2 

10 

2 

26* 

2 

8.4 

2 Jan 1970 

1 

3* 

1 

1 

2 

6.6 

27 Dec 1970 

1 

3* 

1 

1 

2 

1.6 

18 Jan 1971 

3 

2 

68* 

' 6 

19 

I9i6* 

Average 

3.5 

5.0 

8.6* 

6.'2 

4.2 

5.5 


* Relative maximum 


TABLE 18 

Solar Flare In5>ortance (s 2) Vs. Elapsed Days 
Prom Start of Warming 


Elapsed 

Days 

Start Date 



-3 


-4 -5 


3 Jul 1969 2 

I Aug 1969 

II Aug 1969 2 

17 Dec 1969 

17 Dec 1970 



2 



TABLE 


19 

Solar Flare Importance 2) Vs. Elapsed Days 
From Date of Major Warming Peaks 


^ Elapsed 
Days 

Peak Date 



-2 


-3 


-4 -5 


30 Apr 1969 

5 May I969 
8 Jul 1969 
11 Jul 1969 
13 Jul 1969 

6 Aug 1969 

8 Aug 1969 

16 Aug 1969 

26 Dec 1969 

29 Dec 1969 
4 Jan 1961 

8 Jan 1971 

11 Jan 1971 


2 

2 


2 

2 


2 

2 


2 2 

2 3 


2 



2 

2 2 
2 2 
2 2 


2 

2 


t 
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TABLE 
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Solar Flare Importance 2 ) Vs, Elapsed Days 
From Date of Minor Warming Peaks 


Elapsed 

Days 


Peak Date 

0 

-1 

-2 

-3 

-4 

-5 

19 Apr 1969 

2 

2 





23 Apr 1969 

2 


3 

2 

2 

2 

12 May 1969 

3 






3 Jul 1969 

2 

. 

. 2 




18 Jul 1969 

2 






22 Jul 1969 





2 


20 Aug 1969 







13 Dec 1969 







18 Dec 1969 


2 



3 


20 Dec 1969 


2 


2 



2 Jan 1970 

3 






27 Dec 1970 







18 Jan 1971 
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TABLE 
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A Index Vs. Elapsed Cays From Start of Warmings 


Elapsed 

Days 

Start Date 


3 Jul 1969 3 6 

I Aug 1969 3 6 

II Aug 1969 4 6 

17 Dec 1969 4 11* 

17 Dec 1970 3 4 

Average 3.4 6.6 


* Relative maximum 


-2 

-3 

-4 

Average 

17-^ 

5 

4 

7.0 

13* 

2 

4 

5.6 

11* 

9 

8 

7.6 

5 


3. 

5.4, 

14 

65* 

4 

18.0* 

12.0 

17.0* 


8.7 
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TABLE 22 

Ap Vs. Elapsed Days From Date of Major Warming Peaks 


Elapsed 

Days 


Peak Date 

0 

-1 

-2 

.3 

-4 

Average 

30 Apr 1969 

27 

12 

60* 

11 

6 

23 . 2 * 

5 May 1969 

10 

8 

15 

28 * 

8 

13.8 

8 Jul 1969 

5 

6* 

4 

2 

3 

4.0 

11 Jul 1969 

6 

7 * 

7 * 

5 

6 

6.2 

13 Jul 1969 ■ 

V \-*~ 

12 

6 

7 

7 

9.2 

6 Aug. 1969 

5 

7 

13 

14 * 

5 

8.8 

8 Aug 1969 

9 

8 

5 ' 

7 

13 * 

8.4 

16 Aug 1969 

4 

4 

6 

7 

21* 

8.4 

26 Aug 1969 

15 

4 ' 

6 

11* 

5 

8.2 

29 Aug 1969 

4 

, 6 

21* 

15 

4 

9.8 

■4. Jan 1971 

16 

33 * 

19 

6 

2 

9.2 

8 Jan 1971 

0 

1 

4 

9 

16* 

6.0 

11 Jan 1971 

8* 

6 

1 

0 

1 

3.2 

Average 

9.5 

8.8 

12.8* 

9.^ 

7.3 

9.1 


* Relative maximum 




TABLE 23 

Ap Vs. Elapsed Days From Date of Minor Warming Peaks 


Elapsed 


Days 


Peak Date 

t 

0 

-1 

-2 

-3 

-4 

Average 

19 Apr 1969 

5 

15 

21* 

17 

14 

14 , 4 * 

23 Apr 1969 

5 

^ 10* 

6 

8 

5 

6.8 

12 May I969 

8 

5 

10 

11* 

5 

7.8 

3 jui 1969 

3 

6 

17* 

5 

4 

7.0 

18 Jul 1969 

■ 3 ■' 

4 

10 

5 

13* 

7.0 

22 Jul 1969 

7 * 

. 4 

3 

2 

3 

3.8 

20 Aug 1969 

7 

13 * 

7 

6 

4 ' 

' 7.4 

13 Dec 1969 

3 

5 

11 

8 

12* 

7.8 

18 Dec 1969 

- 3 

4 

11* 

5 

4 

3.4 

20 Dec 1969 

3 

4 

3 

4 

11* 

5.0 

2 Jan 1970 

30 * 

8 

3 

2 

4 

9.4 

27 Dec 1970 

8 

4 

4 

14 * 

7 

7., 4 

18 Jan 1971 

15* 

6 

7 

8 

9 

9io 

Average 

7.7 

6.8 

8.7* 

7.3 

5.9 

7.3 
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TABLE 24 


Index Vs. Elapsed Days From Start of Warming 


Elapsed 

Days 


start Date 

0 

-L . 

-2 

-3 . 

-4 

Average 

3 Jul 1969 

12 


7 

18 

18 

11.8 

1 Aug 1969 

6 

-1 

12 

-1 

-18* 

—0*4 

11 Aug 1969 

15 

-6* 

-1 

4 

16 

5.6 

17 Dec 1969 

0 

-6* 

10 

13 

8 

5.0 

17' Dec 1970 

-22 

-35 

-56 

-68* 

10 ' ■ 

-34. 2"^ 

Average 

+2.2 

-8.8^ 


-6.8 ■ 

6.8 

-2.4 


* Relative maximum 
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TABLE 25 

^st Date of Major Warming PeeOcs 

Elapsed 


Days 


Peak Date 

0 

-1 

-2 

-3 

-4 

Average 

30 Apr 1969 

-39* 

-31 

-38 

5 

4 

-19.8* 

5 May I969 

-4 

-9 

-17 

-12 

-21* 

-12.6 

8 Jul 1969 

19* 

24 

29 

24 

19* 

23,0 

11 Jul 1969 

14 

5 * 

8 

19 

24 

14,0 

13 Jul 1969 

2 ^ 

' 7 

14 

5 

8 

7.2 

6 Aug 1969 

4 

-4 

- 13 * 

9 

21 

3.4 

8 Aug 1969 

4 

16 

4 

-4 

- 13 * 

1.^ 

16 Aug 1969 

9 ' 

9 

0 

-4* 

-3 

2,2 

26 Aug 1969 

18 

9 

11 


12 

10.0 

29 Aug 1969 

6 

-12 

-22* 

18 

9 

0.0 

4 Jan 1971 

-8* 

-8* 

7 

5 

-1 

-0.8 

8 Jan 1971 

9 

3 

-2 

-3 

-8* 

-0.2 

11 Jan 1971 

8 

19 

IT 

9 

3 * 

9.6 

Average 

3,2 

2.2 

-0.2* 

5.8 

4.2 

3.0 


* Relative maximum 
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TABtE 
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D Vs. 
st 

Elapsed Elays From Date 

of Minor- Warming 

Peaks 


Elapsed 

Date 

Peak Date 

0 

-1 

-2 

-3 

-4 

Average 

19 Apr 1969 

-10 

-27 

-33^ 

-18 

-22 

-22.0* 

23 Apr' 1969 


15 

-1 

-6 

-10* 

-0.4 

12 May 1969 

14 

3* 

3* 

3* 

10 

6.6 

3 Jul 1969 

la 

4* 

7 

18 

18 

11.8 

18 Jul 1969 

10 . - 

6 

6 . , „ 

. 6 .. 

0* 


22 Jul 1969 

6* 

25 

25 

12 

10 

15.6 

20 Aug 1969 

- 5 * 

2 

11 

-3 

9 

2.8 

13 Dec 1969 

8 

-1 

- 5 * 

2 

9 

2.6 

18 Dec 1969 

10 

0 

-6* 

10 

13 

5.4 

20 Dec 1969 

20 

14 

10 

0 

-6* 

7.6 

2 Jan 1970 

-16* 

9 

13 

8 

6 

4.0 

27 Dec 1970 

-a 

-5 

-8 

-21* 

-2 

-7.6 

18 Jan 1971 

11 

12 

4 

- 5 * 

1 

4.6 

Average 

4.8 

4.4 

2.0 

0.5* 

2.7 

2.9 


* Relative loaximuin 
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TABLE 27 


He II Index Vs. Elapsed From'Start of Warming 


Elapsed 


Days 


Start Date 

0 

-1 

-2 

-3 

-4 

Average 

3 Jul 1969 

82 ^ 

77 

77 

77 

77 

78 

1 Aug 1969 

80 * 

70 

71 

71 

75 

73 

11 Aug 1969 

90 

90 

90 

90 

90 

90* 

Average 

84 * 

79 

79 

79 

81 

80 


* Relative maximum 
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TABLE 


28 


He II Index Vs. Elapsed Days From Date of Major Warming Peaks 


Elapsed 


Days 


Peak Date 

0 

-1 

-2 

-3 

-4 

Average 

30 Apr 1969 

83 

86 

. 88* 

87 

86 

86 

5 May 1969 

85 

84 

84 

85 

86* 

85 

8 Jul 1969 

88 

83 

86 

90* 

90* 

87 

11 JvO . 1969 

87 

89 

91* 

88 

83 

88 

13 Jul 1969 

89 . - 

89 

87 . . 

89 

91* 

89* 

6 Aug 1969 

88 

88 

93* 

86 

83 

88 

8 Aug 1969 

90 

88 

88 

88 

93* 

89*^ 

16 Aug 1969 

75 

77 

82 

85 

87* 

81 

26 Aug 1969 

85 


83 

TO 

70 

79 

29 Aug 1969 

95 ^ 

86 

86 

86 

85 

88 

Average 

86 

86 

87* 

. 85 

84 

86 


* Relative maximum 
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TABLE 29 


He II Index Vs. Elapsed Days From Date of Minor Warming PeaJts 


Elapsed 

Days 


Peak Date 

0 

-1 

19 Apr 1969 

90 

90 

23 Apr 1969 

85 

89 

12 May 19 ^ 

93* 

90 

3 Jul 1969 

90 ^^ 

86 

18 Jul 1969 

84 

84 

22 Jul 1969 

77 

81 

20 Aug 1969 

69 

71 

Average 

■84 

84 


* Relative maximum 


-2 

-3 

-4 

Average 

91 

91 

92* 

91* 

92 if 

90 

90 

89 

88 

88 

85 

89 

82 

77 

77 

82 

85* 

85* 

85* 

85 

82 

82 

84 * 

81 

72 

73 , 

75 * 

72 

84 

84 

84 

84 
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TABLE 30 


SID Index Vs. Elapsed Days from Start of Warming 

j 

Elapsed 

Days 

Start Date 0 ,-1-2 -3 -4 

3 Jul 1969 3 2-^ 

1 Aug 1969 2 3 2^ 

11. Aug 1969 3" 

17 Dec 1969 2 2 

17 Dec 1970 
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TABLE 31 


SID Index Vs. Elapsed Days From Date of Majdr Warming Peaks 


Elapsed 

Days 

Peak Date 0 

30 Apr 1969 

5 May 1969 2 

8 Jul 1969 

11 Jul 1969 
13 Jul 1969 

6 Aug 1969 
8 Aug 1969 
16 Aug 1969 
26 Aug 1969 
29 Aug 1969 
4 Jan 1971 
8 Jan 1971 
11 Jan 1971 
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TABLE 32 


SID Index Vs, 

Elapsed 

Days 

Peak Date 

Elapsed Days From Date of Minor Warming Peaks 

0 -1 -2 -3 -4 

19 Apr 1969 

2" 2"^ 

2+ 

23 Apr 1969 

2^ 

2*^ 2" 

12 May 1969 



3 Jul 1969 

3 2*^ 


18-jta 1969 


... — ; ' ' ■ ■ 

22 Jul 1969 ■ 



20 Aug 1969 

2 


13 Dec 1969 

2 3 " 


18 Dec 1969 

2 

2 

20 Dec 1969 


2 

2 Jan 1970 


2 

27 Dec 1970 

2 


18 Jan 1971 

2^ 

2 




TABLE 33 


Magnetic Condition Index Vs. Elapsed Days From Start 

of Warmings 

Elapsed 

Days 

Start Date 

0 -1 -2 

-3 -4 

3 Jul 1969 

S 


1 Aug 1969 

S 


11 Aug 1969 

s 

S 

17 Dec 1969 



17 Dec 1970 


.S . 
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TABLE 34 

Magnetic Condition Index Vs. Elapsed Days frcam Date of Major Warming Peaks 


Elapsed 

Days 

Peak Date 0 -1 -2 -3 


30 Apr 1969 

^ S s 


5 May 1969 

S 


8 Jul 1969 
11 July 1969 
13 July 1969 

S 


6 Aug 1969 

D 


8 Aug 1969 

s 


16 Aug 1969 


S 

26 Aug 1969 

3 


29 Aug 1969 

D S 


4 Jan 1971 
8 Jan 1971 
11 Jan 1971 




D ■ Disturbed Day 
S ■ PrincipeG. Magnetic Storm 
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TABLE 35 


Magnetic Condition Vs. Elapsed Days From Date of Minor Warming PeaXs 


Elapsed 

Iteite 

PeaX Date 


0 




-4 


19 Apr 1969 
23 Apr 1969 

12 May I969 
3 Jul 1969 
18' Jul 1969 

22 Jul 1969 

20 Aug 1969 

13 Dec 1969 
18 Dec 1969 
20 Dec 1969 
2 Jan 1970 
27 Dec 1970 
18 Jan 1971 
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TABLE 36 

Average Value, of Indices Vs, Warming Days, Wohwaiming Days, Total Days, and 
For 0-4 Days Freceding S^art of Warming, Major Warming Peaks, and Minor 

Warming Peaks 


Index 

Warming 

Days 

Non- 
Warmlng 
Bay s’ , 

Total- 

Days 

Start of 
Warming 

Major 

Warming 

Peaks 

Minor 

Warming 

Peaks 

AE 

170 (H)* 

178(H) 

175 

141(B) 

198 (A) 

149 (B) 

A 

P 

9 . 0 (B) 

10.8(A) 

10.0 

8 . 7 (B) 

9 . 1 (B) 

7.3(B) 

®st 

+0.7(B) 

-2.2(A) 

- 0.6 

-2.4(A) 

3 . 0 (B) 

2 . 9 (B) 


* N = Normal value 

B = Below-normal value 
A B Above-normal value 
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Fig. 1 March of Maximum Zonal Radiance, SIRS 
Channel 8, April 15-May 14, 1969 
(Circled points are hemispheric maxima) 
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APRIL 1969 MAY 1969 


Fig. 2 March of Maximum Radiance, SIRS A Channel 8 
April 15-May 14, 1969 



RADIANCE(ERGS CM"^ S"' STER"’ CM" 



JULY 1969 


S," ® •'"ir 1969 

Solid curve is maximum radiance 

temperatures from 

qmf+K' and Mclnturff, 196Q; 

Smith, et al, 1970) » j 
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Fig, 4 March of Maximum SIRS A Channel 8 Radiance, Jtay 25-AugustJ24, 1969 
Southern Hemisphere i 

Solid curve is maximum radiance. 

Dashed curves are 10- and 30-mb temperatures from Mirny, Antarctica 
(Finger and Mclnturff, 1969, Smith, et al, 1970) 
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Fig. 5. Change in radiances(ergs/) m:i®U3.3 SIRS(669.3, cm"^). from 
25 June 1969 to 10 July 1969.| (Fritz and Sox^ee, 1970) 
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Fig. 6, Vfoighting functions for SIRS channels tWark, 1970) 
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Pig. 7 30-mb analysis for 17 August, 1969. Units: geo- 

; potential meters and degrees Celsius. (Miller; et al, 1970) 



KEPBODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 



Figured. 2 mb analysia for 24 December 1969. Filled circle represents an observation for the analysis day, 
an open square one t^t is up to two days from analysis day. A barb on the wind vector represents 10 kt, a 
filled pennant 50 kt. Flight is presented in geopotentiai decameters, temperature in °C. 

"Mi l ler , «t a l, 1 972 ) 
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JANUARY 14, 1970 


FigurcXX Same as Fig 6 for 14 January 1970. 
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SIRS CHANNEL 8 RADIANCE 
er«6 cro“^sec"^ster*-^/cm' 

* JAN. 5, 1970 , 

Fig ^-4-. Channel 8 radiance map for 5 January 1970. Principal hlgh-radtance center (66 ergs/) is located over NE Greenland coast. 
A secondary high (61.5 ergs/) was generated from main system in polar area, travelling westward to Far East 
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Fig. 15 Northern Hemisphere 10-ii^ km) Chart for Becember 26, 1969 

Maximum warming center (266^) is located at llO^E, 57*^N. 

V\ ’’ / ’^, ■ V~i > — i \. •••*-' j\ • m' \ »»f u.ij o,. 


i, ^ 








Fig. 16 Movement of Warming Center at 10 mb from 25 to 30 December 1969 

Date is beneath locating cross. 

Circles locate warming center at 2 mb (43 km) 



Fig. 17 March of Maximum Daily 10-mb Temperature at Latitudes Above 40"^5■ 
Between December 5, 1969 and January 14/1970 

Dashed curve indicates march of 2-mb (43 km) Temperature 


2 MB TBMBEHATURE (°K) 




\r~/ /*\ 








Fig. a, Teraparatures at 30 mb(solid lines) and radiances of SIRS channel ^(dashed 
lines) for 7 Jan* 1971. b, Temperatures at 45 km(heavy lines) and radiances(thin 
lines) of channel A of SCR for 6 Jan. 1971, based on all rocket observations available 
within 3 days before and after map date, c, temperatures at 60 km for 6 Jan, 1971 
(with the notation of Fig. b)* (Labitzke, 1972) 
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20 25 30 31 5 10 15 20 24 

DECEMBER 1970 JANUARY 1971 


Fig. 20 March of Maximum Daily 10-mb Temperature Above dO^N Latitude 
M From December 15 » 1970 to January 24 f 1971. 

I ■ . ’ 

o Dashed curve indicates the march of 2-mb (45 1cm) high-latitude 

^ zonal average temperature. 




Fig. 22 Temperature-Height Brofile in the Vicinity of a 
Varmtng, Derived from Constant Pressure Charts 
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Fig. 23 


Locations of Warm and Cold Centers from 25 to 60 km 
on January 6, 1971 in the ITorthem Hemisphere 

^ kilometers is indicated below locations. 
Triangles indicate positions of warm centers at 31 km 
Square indicates position of low. 









Fig. 25 Locations of Temperature Change Centers from 25 to 60 km 
January 6, 1971- 

Temperature increase centers indicated by crosses. 
Temperature decrease centers indicated by circles. 
Altitude in kilometers shown below circle or cross. 
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Pig. 26 Maximum Temperature Increase at Mtitudes Between 
25 and 60 Kilometers at 78°N, 70^ During 1970-71 
Warming . 
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AN HyPOTHESIS FOR THE RELATION BETWEEN 
AURORAL ELECTROJEPS AND STEUiTOSPHERIC WARMINGS 

Albert D. Anderson 


' The following explanation is presented to show how warming events 
at' 10 mb may be related to magnetic activity, especially AE, with a time 
delay of about two days. First, the properties of auroi*al electrojets, 
the winter circumpolar vortex and planetary waves will be discussed. 

The Auroral Electrojets (AEJ) 

Thermospheric heating is definitely , correlated with magnetic . 

activity as indicated by A^ and Dst (Anderson, 19?3). Cole (197I) indicates 
that the heating mechanism is Joule dissipation associated with auroral 
electrojets (AEJ) that results in a- global input of heat into the lower 
thermosphere during strong geomagnetic disturbances that is larger than that 
due to solar extreme ultraviolet radiation. The marked correlation 
exhibited between the magnetic indices, especially the AE index, and strato- 
spheric warmings suggests that AEJ heating strongly influences the lower 
atmosphere also, with a time delay of about 2 days. 

AEJ are intense electric currents, associated with polar magnetic 
disturbances, that flow in the auroral oval in the lower thermosphere 
(Akasofu et al . , 19^5) . The westward AEJ grows to approximately 10^ a.mperes 
during a magnetic substorm due to an increase in the conductivity of the 

E-layer of the ionosphere in the midnight to dawn sector of the auroral oval. 

3 

The eastward flowing AEJ has a strength of about 10 amperes in the late 
evening sector. Cole (I962) has estimated that 5 ergs cm’^s"^ could be 
deposited in the auroral zone by means of Joule heating associated with the 
AEJ at moderately disturbed times and an order of magnitude more during 
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large geomagnetic storms. The height of majcimum joule heating lies between 
100-150 Ion (Blumen and Hendl, I969). Cole (l9Tla) indicates that major 
AEJ heating of the thermosphere during strong geomagnetic disturbances could 
perturb the atmosphere down to 100 km. 

Winter Circumpolar Vortex 

It is well known that the winter circulation in the troposphere is 
dominated by a cold low centered near the pole that becomes more pronounced 
at higher altitudes and is surrounded by a band of strong westerly winds 
(Labitzke, I972). Evidently, this vortex extends up through the stratosphere 
and mesosphere to the base of the thermosphere near 100 km. The normal wes- 
tefly zonal flow in the meteor region (75-IO5 km) is consistent with a 
large-scale cyclonic system centered near- the pole (Hook-, 1972). However, 
during warmings there, is evidence Indicating substantial circulation and 
temperature changes In the mesosphere (Quiroz, I969; Hook, I970). During 
the midwinter stratospheric warmings for 1967-68 and I970-7I, both the 
meridional and zonal winds at College, Alaska at 75-IOO km reversed (Hook, 
1972). The reversal of the zonal flow during the warmings is consistent 
with an antlcyclonic system. During the warming of I970-7I, the . prevailing 
flow reversed from about 3l December 1970 to- 7 January I971. Hook (I972) 
concludes that the warming in the upper mesosphere started approximately 10 
days before the warming at the 10 mb level. From this, he states, that the 
circulation breakdown during this warming affected the upper mesosphere 
before it affected the mid-stratospheric level. 

Plane-tary Waves and the Critical Level (CL) 

Planetary-scale disturbances can propagate upward from the tropo- 
sphere into the higher atmosphere if the prevailing winds are moderate 
wes^^erlies (Charney and Drazin, I961). Planetary waves occur in the strato- 
sphere (Finger et al . , I966), reaching as high as the thermosphere (Newell 
and- Dickinson, I967). The waves cannot propagate upward in easterly winds 
so that when "the waves reach a critical level (CL), the level at which the 
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zonal velocity is zero, the wave amplitude decreases rapidly and the wave 
energy flow vanishes (Dickinson, I968), If the level is horizontal, 
meridional heat transfer by the waves is constant in the westerly layer 
and jumps to zero at the CL. Hence, the temperature disturbance amplitude 
is very large near the CL, The easterly winds near the CL are accelerated 
by the waves, thereby lowering the level (Matsuno, I971). Thus, the warm- 
ing and wind reversal shift downward. The easterly layer spreads downward 
first very rapidly, then slowly, because the air density inpreases as the 
layer lowers. The warmest layer moves down with the descent of the CL. 

Since the waves transport heat from lower to higher latitudes, the tempera- 
tures at the hi^ latitudes is always higher than at the lower latitudes, 
except initially. Intense warming of the polar region is expected to occur 
just below the CL, while the temperature may fall in the low-latitude region. 

An Hypothesis for Thermospheric-Stratospheric Interaction (Warmings) 

Taking into account the information concerning the AEJ, Cold Low, 
Pl^etary Waves and CL, an attempt will now be made to outline the dynamics 
of a stratospheric warming; 

Strong AEJ heating occurring in the auroral zone hear 100 km is 
conducted down into the mesosphere, destroying the cold low from above and 
replacing it with a warm anticyclone. This results in the formation of the 
CL in the upper mesosphere. Johnson and Gottlieb (I970) have shown the 
important role played by the eddy transfer of heat in the 60-120 km region; 
at higher altitudes heat transfer between levels is dominated by molecular 
conduction, and at lower levels by radiation exchange. Downward motion in 
the upper mesosphere will be accompanied by compressional heating, thus 
enhancing the AEJ heating. Normally, a strong planetary wave propagated 
upward from the troposphere would be dissipated in the thermosphere. However, 
once the CL is formed, the wave interacts with the easterly winds at the CL 
in the upper mesosphere, heats the air there, enhances the CL and lowers it. 
Once this process is started, the lowering of the CL will continue, as long 
as the strong wave exists. If the wave is strong enough and persists long 


III-3 



enough, the CL will lower down into the stratosphere, even though the AEJ 
is no longer heating the upper mesosphere and cooling is taking place there 
(Labitzke, 1972). When the critical level descends to near 4-5 km, "explosive" 
Warming starts, probably because the wave energy density attains a maximum 
in the middle and low stratosphere (Matsuno, 1970)- The length of time it 
takes for a strong CL to descend from the upper atmosphere into the strato- 
sphere is unknown. It could take arywtiere from a day to a week or more. 

Apparently, AEJ heating is not continuous during warmings. Evidently, 
there Is no significant difference between the average AE index for warm- 
ing days and nonwarming days- The AE index varies greatly from hour to 
hour. Therefore, in order for the CL to form, the strongly fluctuating AEJ 
heating must proceed in such a fashion as to progressively destroy from 
above the cold low in the upper mesosphere. Once the' CL is established,"the" 
interplay among the amount and duration of AEJ heating and amplitude of waves 
from the troposphere determines its Initial strength and altitude, and this 
in turn determines the starting altitude and initial extent of the warming. 


Although the existence of a strong CL descending in the mesosphere 
can explain some features of the warmings occurring there, a strong CL 
cannot explain most of the characteristics of the warmings in the stratosphere. 
For example, the warming in I969-I970 started about December I7, I969 and 
ended about January 5, 1970. Hence, it affected the 10 mb altitude for 
seWral weeks, a time interval typical of major warmings. Yet the start of 
warmings and many of the major and minor peaks appear at 10 mb many days 
before the "explosive" warming in the stratosphere. Apparently, pulsations 
of warming in the stratosphere precede a major warming event (Hook, 1972). 
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During the December I969 - January I970 and Deceiriber I970 - January I97I 
Northern-Heminphere warmings many small eastw^d-moving warm (closed Isothenn) 
centers appeared before the major warmings at midlatitudes ( 30 *^W to 50°N) 
on the periphery of the low circulation. Evidently, the warming builds up 
in the stratosphere, at least in part, by successive heat pulses. Although 
heat pulses are usually observed before a major warming, they do not appear 
to constitute a necessary condition for the development of a major warming 
(Quiroz, 1973). The generation of these pulses appear to start at high 
mesospheric levels in early winter (Theon and Smith, I970). It ;has already 
been noted that the major fluctuations in the AE index indicates that the 
AEJ heating is not a continous process. Also, the fact that there is roughly 
the same time delay (about two days) between the maximum AE index and start 
of warmings, and major and minor peaks at 10 mb indicates that the 
AEJ, heating near 100 km could be responsible for the warming phenomena at 
10 mb. ■ 

The' following tentative expleoiatlon is offered for the pulsations of 
warming at 10 mb: It is assumed that (a) a weak CL exis-ts in the upper 
mesosphere only when AEJ heating Is above a certain intensity, and (b) waves 
are reflected (Hines, 1972 ) "by the weak CL down into the stratosphere where 
they interact with upward propagating waves to release heat in the region 
where the wave energy density attains a maximum. A weak CL implies the 
absence of very strong waves, although weak waves may be present. Since the 
AEJ heating is fluctuating, the weak CL probably would have a sporadic exis- 
tence in the upper mesosphere. Thus, even though the waves are persistent, 
the reflection of waves from the upper mesosphere fluctuates directly with 
AEJ heating. The reflected waves take roughly two days to descend to 10 mb. 

N 

Some of the heating characteristics that appear at 10 mb may depend 
also on the horizontal shape of the CL. Initi ally , the CL should form in 
the auroral zone. However, the shape of the CL probably is not symmetric 
with the auroral zone. The initial CL shape will depend on the interaction 
between the newly-developing anticyclonic circulation resulting from CL 
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heating with the strong western^ winds around the cold polar low. Some 
of the heat created at the CL prohably will not remain there; a portion' 
be transported vertically and some may be ducted horizontally 
(Friedman, I966; Cowling et al . , I971), 

The horizontal extent of the CL in the upper mesosphere cannot be 
determined at present because of lack of data. Perhaps, some idea of 
its distribution in the stratosphere can be obtained by examining available 
Northern Hemisphere stratospheric Isotherm patterns on a date when a major 
warming appeared there. Figure 1 displays a disturbed temperature map 
for the Northern Hemisphere at 50 km, January 6 , I971, a major warming 
date. The level of maximum warming was between’ 4-0 and 
50 km on this date. Gwling took place above the warming area (Labltzke, 
1972). Hence, Figure 1 may represent heating near the altitude of the CL. 
The heating pattern In Figure 1 appears quite similar to the theoretical 
disturbed geopotential height distribution at 45 km computed by Matsuno 
(1970) for^ the January I967 warming. Yet there Is one marked difference. 

In Figure 1 the zone of maximum warming closely parallels the auroral 
zone for' over half of its perimeter on the east. This is not true for 
■the disturbed geopotential height distribution. During a geomagnetic 
storm, the auroral zone would expand to approach even more closely to 
the zone of maximum heating in Figure 1 . 

The close association of the zone of maximum heating with the 
auroral zone in Figure 1 indicates that at least part of the original 
heating pattern established at the CL in the upper mesosphere has per- 
sisted down to 50 km. Willett (I968) has pointed out the existence of a 
circumpolar zone of maximum warmth in the arctic winter stratosphere, 
centered on the geomagnetic pole. All of the observed stratospheric 
warmings from 1951 to I966 developed in the equatorward dip of the 
auroral zone. In Figure 1 , the heating extends to midlatitudes at longi- 
tudes near the equatorward dip of the auroral zone; at 45 °N the zone 
of maximum heating is at the same longitude ( 290 °E) as the dip. It was 
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pointed out that many small warming centers appear initially at 10 mb at 
midlatitudes (30 W to 50 °N) and move east on the periphery of the lov 
circulation. Also, it is observed that many of the centers initially 
appear at longitudes east of the auroral trough (290°E) in the 300°E to 
360 E sector. Although the zone of maximum warming in Figure 1 extends 
from high to mid-latitudes ^ the heat at midlatitude may be noticeable 
because it appears as local temperature maxima 'at 10 mb. Wo small, 
closed warm isotherm patterns develop at higher latitudes because the 
heat is mixed with the cold ascending currents of the low. This explana- 
tion is based on the assumption that the heat pattern in Figure 1 is 
representative of a weak CL pattern in the upper mesosphere. Consequently, 
the heat at mid-latitudes results from the Interaction of the reflected 
wave with the upward wave in the mid-stratosphere. 

Implications for Other Atmospheric Phenomena 

The magnetic Indices are correlated somewhat with each other, A^ 
contains contributions from at least two major sources; the auroral electro 
Jet (AE) and the ring current (Dst). Therefore, AEJ heating may be the 
key to explaining other puzzling correlations between magnetic activity 
and lower atmosphere behavior. Studies (Levine et al . . 1974) have indi- 
cated that meteorological responses in the lower atmosphere occur within 
2 or 3 days after geomagnetic activity and that these responses are most 
pronounced in winter. If AEJ heating affects the circulation in the meso- 
sphere and stratosphere, undoubtedly it will influence conditions in the 
troposphere too. In other words, during a major warming, the whole 
circulation in the winter atmosphere, extending from the troposphere to 
the thermosphere, probably is affected. If AEJ activity affects the lower 
atmosphere, then, inasmuch as AEJ heating is controlled by solar (plasma) 
activity, it follows that variations in solar activity affect the 
weather. Consequently, a new link between solar activity and the lower 
atmosphere may have been uncovered. 
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